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ABSTRACT 


The on-off reaction jet control system is often used for attitude and 
orbital maneuvering of various spacecraft. Future space vehicles such as the 
orbital transfer vehicles, orbital maneuvering vehicles, and space station will 
extensively use reaction jets for orbital maneuvering and attitude stabilization. 
The proposed robust fuel- and time-optimal control algorithm is used for a 
three-mass spring model of flexible spacecraft. A fuel-efficient on-off control 
logic is developed for robust rest-to-rest maneuver of a flexible vehicle with 
minimum excitation of structural modes. 

The first part of this report is concerned with the problem of select- 
ing a proper pair of jets for practical trade-offs among the maneuvering time, 
fuel consumption, structural mode excitation, and performance robustness. A 
time-optimal control problem subject to parameter robustness constraints is 
formulated and solved. 

The second part of this report deals with obtaining parameter insensi- 
tive fuel- and time-optimal control inputs by solving a constrained optimization 
problem subject to robustness constraints. It is shown that sensitivity to mod- 
eling errors can be significantly reduced by the proposed, robustified open-loop 
control approach. 

The final part of this report deals with sliding mode control design for 
uncertain flexible structures. The benchmark problem of a flexible structure 
is used as an example for the feedback sliding mode controller design wdth 
bounded control inputs and robustness to parameter variations is investigated. 







ACKNOWLEDGEMENTS 


This research was supported by the NASA Johnson Space Center 
through the RICIS program of the University of Houston at Clear Lake. The 
authors would like to express special thanks to Dr. Kenneth Cox and Dr. John 
Sunkel of the NASA JSC for sponsoring this research. 







TABLE OF CONTENTS 


LIST OF TABLES 

LIST OF FIGURES 

Chapter 

1 . Introduction 

2 . Robust Time-Optimal Control of Uncertain Flexible Spacecraft 

2. 1 Introduction 

2.2 Time-Optimal Rest-to- Rest Maneuver Control 

2.2.1 Problem Formulation 

2.2.2 Case 2-1 with a Two-Sided Control Input 

2.2.3 Case 2-2 with Two One-Sided Control Inputs 

2.2.4 Case 2-3 with Two One-Sided Control Inputs 

2.3 Robust Time-Optimal Control 

2.3.1 Problem Formulation 

2.3.2 Case 2-1 with a Two-Sided Control Input 

2.3.3 Case 2-2 with Two One-Sided Control Inputs 

2.3.4 Case 2-3 with Two One-Sided Control Inputs 

2.4 Summary 

2.5 Conclusions 

3 . Robust Fuel- and Time-Optimal Control of Uncertain Flexible Spacecraft . 

3 . 1 Introduction 

3.2 Fuel- and Time-Optimal Control 

3.2.1 Problem Formulation 

3.2.2 Rigid Body Fuel- and Time Optimal Control 


Page 

vii 

viii 

1 

5 

5 

7 

7 

12 

14 

16 

18 

18 

19 

21 

23 

25 

26 
36 
36 
38 

38 

39 








Chapter 


Page 

3.2.3 Parameter Optimization Problem 40 

3.2.4 Case 3-1 42 

3.2.5 Case 3- lb 43 

3.2.6 Case 3-2 44 

3.2.7 Case 2-1 45 

3.2.8 Case 2-lb 47 

3.3 Robust Fuel- and Time-Optimal Control 47 

3.3.1 Problem Formulation 47 

3.3.2 Case 3-1 49 

3.3.3 Case 3-2 51 

3.3.4 Case 2-1 52 

3.4 Conclusions 54 

4. Sliding Mode Control Design for Uncertain Flexible Structures 70 

4.1 Introduction 70 

4.2 Variable Structure System 71 

4.3 Sliding Mode Control Theory 72 

4.4 Model of an Uncertain Flexible Structure 88 

4.5 Sliding Mode Controller Design 89 

4.6 Conclusions and Issues 92 

5. Conclusion 101 

REFERENCES 104 

APPENDIX 

A. Numerical Solution of Nonlinear Constrained Optimization Problem 109 


vi 






UST OF TABLES 


Table 

1 


Summary of the results 


Page 

26 


vii 




LIST OF FIGURES 


Figure 

2.1 Three-mass-spring dynamical system 6 

2.2 Pulse sequences ^ 

2.3 Responses of time-optimal control Case 2-1 28 

2.4 Responses of time-optimal control Case 2-2 29 

2.5 Responses of time-optimal control Case 2-3 30 

2.6 Responses of robust time-optimal control Case 2-1 31 

2.7 Responses of robust time-optimal control Case 2-2 32 

2.8 Responses of robust time-optimal control Case 2-3 33 

2.9 Responses of robust time-optimal control Case 2-3 34 

2. 10 Responses of robust time-optimal control Case 2-3 35 

3.1 Case 3-1 ^ 

3.2 Case 3-2 56 

3.3 Case 2-1 57 

3.4 Pulse sequences ^7 

3.5 Responses to fuel and time-optimal control input (Case 3-1) 58 

3.6 Responses to fuel and time-optimal control input (Case 3- lb) 59 

3.7 Responses to fuel and time-optimal control input (Case 3-2) 60 

3.8 Responses to fuel and time-optimal control input (Case 2-1) 61 

3.9 Responses to fuel and time-optimal control input (Rigid Body Solution) . 62 

3. 10 Responses to fuel and time-optimal control input (Case 2- lb) 63 

3. 1 1 Responses to robust fuel and time-optimal control input (Case 3- 1 ) 64 

3. 1 2 Responses to robust fuel and time-optimal control input (Case 3-2) 65 

3.13 Responses to robust fuel and time-optimal control input (Case 2-1) 66 

3. 14 Responses to fuel and time-optimal control input (Rigid Body Solution). . 67 

viii 




Figure Pa g e 

3.15 Responses to robust fuel and time-optimal control input (Case 2-1) 68 

3.16 Responses to robust fuel and time-optimal control input (Case 2-1) 69 

4. 1 Sliding Mode Region 93 

4.2 Sliding Surface 94 

4.3 Math model of two-mass-spring dynamical system 95 

4.4 Response and control of sliding mode controller for the case with 

parameters k = 1, ml = 1, m2 = 1 (nominal case) 96 

4.5 Response and control of sliding mode controller for the case with 

parameters k = 1, ml = 0.5, m2 = 0.5 97 

4.6 Response and control of sliding mode controller for the case with 

parameters k = 1, ml = 0.7, m2 = 0.7 98 

4.7 Response and control of sliding mode controller for the case with 

parameters k = 1, ml = 1.5, m2 = 1.5 99 

4.8 Response and control of sliding mode controller for the case with 

parameters k = 1, ml = 1, m2 = 0.7 100 




IX 





Chapter 1 


Introduction 


In recent years there has been an interest in deployment of large space 
structures in earth’s orbit for scientific and commercial activity. These large 
structures are required to be lightweight and are flexible compared to the bulky 
rigid structures which were deployed earlier. These flexible structures need ac- 
tive control to damp vibrations and to maintain figure. The design of control 
systems for such an application using nonlinear control elements like an on-off 
thruster is a different problem. One problem is controlling these structures for 
rest-to-rest maneuver in minimum time. The well known standard problem of 
bang-bang control has been solved using the minimum principle in [1], The re- 
quirement of time minimization results in bang-bang or pulsed one sided control 
which can be implemented using current on-off actuation technology. Another 
problem would be to minimize a performance index which weighs time against 
fuel consumption for a rest-to-rest maneuver. This is a more difficult problem 
for trajectories originating from the same point which satisfy the first order 
necessary conditions are not unique. Further, the control system should be ro- 
bust to uncertainty in modeling which would occur due to parameter varations. 
The total inertia of the flexible structure is well known but there is uncertainty 
in the mass distribution and stiffness of the flexible structure. These parameter 
variations has been shown to create residual vibration during and/or after ma- 
neuver. The problem is of selecting a proper pair of jets for practical trade-offs 



0 


among maneuvering time, fuel consumption, structural mode excitation, and 
performance robustness is studied. Open-loop study of robust fuel- and time- 
optimal control gives insight into the characterstics of control structures which 
optimizes manuevering time and fuel consumption. Such controls would be 
desirable in closed-loop control system designed to reduce the effect of distur- 
bances and parameter variations. The ultimate goal of the future is to develop 
nonlinear feedback control logic for achieving the robust fuel- and time-optimal 
performance. Most robust linear feedback controllers do not fully utilize the 
available control power in performing a time-optimal maneuver. An attempt 
is made to use the sliding mode control technique to obtain a fast settling time 
and be robust to parameter variations using bounded control inputs. 

The minimum time problem for flexible spacecraft has been addressed 
in several technical papers [2, 3, 4, 5]. In [3], the bang-bang control is shown to 
have an odd number of switches and to be an odd function of time about the 
middle switch. In [3,4,5] statements are made that “in most cases” the time- 
optimal control for a rigid body and n bending modes has 2n + 1 switches. In 
[6] it is proved that when n = 1, there are never more than three switches for 
scalar control. In [7] an optimal open-loop control is obtained for fixed ma- 
neuver time where the intermediate switching time is selected to minimize the 
post-maneuver energy. In [14,15] near-optimal maneuvers to perform the retar- 
getings, which minimize fuel consumption, were also obtained. In minimizing 
settling time, it is shown that time-optimal control is sensitive to parameter 
variations. The residual vibrations occurs as a result of parameter variations. 
In [10] forcing functions are constructed which attenuate residual vibrations. 
In [8] a finite Fourier expansion for the forcing function is used and coefficients 
are picked in such a way so as to reduce the peaks of the frequency spectrum 
at discrete points. This only eliminates a few of the peaks and leaves some 
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modes excited. However this is not time-optimal forcing function. In [9] an 
appropriately shaped torque is derived which not only minimizes residual vi- 
brations but also minimizes the effect of parameter variations which change the 
modal frequencies. These methods use a smooth continuous forcing functions 
having initial and final zero slopes so that the structural modes are not excited. 
But our results show to contrary that the proper sequence of on-off pulse are 
required. 

The sliding mode control theory has been studied in great detail in 
[24]. This nonlinear feedback control is used to achieve accurate tracking in 
presence of disturbance and parameter variations. In [22] continuous control 
laws are incorporated to approximate the idealized sliding mode control. This 
eliminates chattering, the high frequency component of control which occurs 
in idealized sliding mode control. However, these schemes do not use bounded 
control inputs. 

In this thesis the primary objective is to study the feasibility of com- 
puting open-loop on-off pulse control logic for uncertain spacecraft. The robust 
time-optimal control which brings the flexible spacecraft from rest-to-rest ma- 
neuver are obtained. The lumped mass spring model with negligible damping is 
used as a model of a flexible spacecraft. A constrained parameter optimization 
approach is used to solve this problem. The theoretical and practical imple- 
mentation issues inherent to constrained optimization are not elaborated. The 
bang-bang control and one sided control inputs are obtained. The one-sided 
control need not be an odd function about mid maneuver time and it becomes 
a difficult problem lo solve than the standard problem. Many theoretical issues 
(i.e. uniqueness and structure of time-optimal) needs to be resolved. It is shown 
that an accurate mathematical modeling is required for time-optimal solution 
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and that they are sensitive to modeling uncertainty. The parameter optimiza- 
tion problem is formulated and solved with additional robustness constraints. 
Because of properly coordinated pulse sequences the flexible modes are not 
significantly excited during maneuver and the residual response after the ma- 
neuver are well desensitized. The two one-sided controls show “bang-off-bang” 
nature in the time optimal solution. This motivates the problem of computing 
the on-off control logic for fuel- and time-optimal control of a flexible spacecraft 
which is generally a difficult problem to solve. The mass spring model with 
one/two flexible modes are used as examples of flexible spacecraft. The sensi- 
tivity to parameter variations is reduced by using robustness constraints. The 
feedback control method of sliding mode control design for uncertain flexible 
structure for benchmark problem using bounded control inputs is presented. 
The robustness to parameter variations is illustrated for this design. 

The remainder of this thesis is organized as follows. Chapter 2 deals 
with the robust time-optimal control of uncertain spacecraft. In chapter 3, 
the robust fuel- and time-optimal control of uncertain flexible spacecraft is 
presented. Chapter 4 presents the detailed review of the sliding control theory 
and its application to the benchmark problem of flexible spacecraft. Chapter 
5 is the conclusion. Appendix A reviews the numerical solution of nonlinear 
constrained optimization problem. 


Chapter 2 


Robust Time-Optimal Control of Uncertain Flexible 

Spacecraft 


2.1 Introduction 

This chapter deals with the problem of computing open-loop, on-off 
jet firing logic for flexible spacecraft which are sometimes required to maneuver 
as quickly as possible with a minimum of structural vibrations during and/or 
after a maneuver. Most standard time-optimal control approaches to such 
a problem require an accurate mathematical model, and thus the resulting 
solution becomes sensitive to variations in model parameters. 

Expanding on the recent results of [11,12], we further explore the 
robust time-optimal control problem of flexible spacecraft in the face of mod- 
eling uncertainty. In particular, we study the problem of selecting a proper 
pair of jets for practical trade-offs among the maneuvering time, fuel consump- 
tion, structural mode excitation, and performance robustness. A parameter 
optimization approach, with additional constraints for performance robustness 
with respect to modeling uncertainty, is employed to solve such a robust time- 
optimal control problem. However, many theoretical and practical implemen- 
tation issues inherent to constrained parameter optimization problems are not 
elaborated in this chapter. 

A simple math model of flexible spacecraft with a rigid- body mode 
and two flexible modes, as shown in Fig. 2.1, is used to illustrate the concept and 
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methodology. We consider the case in which the structural flexibility and mass 
distribution of the vehicle are quite uncertain, while the total mass (or inertia) 
of the vehicle is well known. Consequently, we focus on the robust control 
problem of flexible spacecraft in the face of modal frequency uncertainty as 
well as mode shape uncertainty. 


u„ u. 


m i m 2 m 


Case #1: hi I £ 1 , u 2 ■ 0 , u 0 
Case #2: OSu^+l.-lSu^u *0 
Case #3: 0 £ Uj £+1, 0, -IS u ^0 


Figure 2.1: Three- mass-spring dynamical system. 

Other robustifled, open-loop approaches, however, attempt to find a 
smooth continuous forcing function (e.g., a versine function) that begins and 
ends with zero slope. The basic idea behind such approaches is that a smooth 
control input without sharp transitions is less likely to excite structural modes 
during maneuvers. On the contrary to such a common notion, the results of this 
chapter indicate that properly modulated, on-off pulse sequences can achieve a 
fast maneuvering time with a minimum of structural vibrations during and/or 




after a maneuver, even in the face of plant modeling uncertainty. 

The remainder of this chapter is organized as follows. Section 2.2 de- 
scribes the standard time-optimal control problem of flexible spacecraft with- 
out modeling uncertainty. A parameter optimization problem is formulated, in 
which the objective function is the maneuvering time. Three cases are explored, 
as illustrated in Fig. 2.1, in order to assess the actuator placement problem for 
time-optimal control of flexible spacecraft with multiple jets. In Section 2.3, 
we investigate the same problem as in Section 2.2 but considering the presence 
of modeling uncertainty. A time-optimal control problem subject to additional 
robustness (or sensitivity) constraints is formulated, and numerical solutions 
for three cases are then compared with solutions obtained in Section 2.2. 

2.2 Time-Optimal Rest-to-Rest Maneuver Control 

2.2.1 Problem Formulation 

Consider a flexible spacecraft described by 

Mx + Kx = Gu (2.1) 

where x is a generalized displacement vector, M a mass matrix, K a stiffness 
matrix, G the control input distribution matrix, and u the control input vector. 

Equation (2.1) is transformed into the modal equations: 

yi + w,yj = i^nUi -f <j> uii2 + ^i3U 3 

$2 + ^>lV2 — 4> 21 U 1 + 4>22 u 2 + ^ 23^3 . . 

Vn + = <t> n \U\ + <An2«2 + 

where yi is the i th modal coordinate, u>, the i th modal frequency, 4> X1 the modal 
input distribution coefficient, and n the number of modes considered in control 



design. Without loss of generality, only three control inputs are considered in 
Eq. (2.2). 

In this chapter, we consider a simple model which is a generic rep- 
resentation of a spacecraft with a rigid-body mode and two flexible modes, as 
shown in Fig. 2.1. Three cases are studied: (i) Case 2-1 with both “positive” 
and “negative” jets placed at body 1, (ii) Case 2-2 with a “positive” jet at 
body 1 and a “negative” jet at body 2, and (iii) Case 2-3 with a “positive” jet 
at body 1 and a “negative” jet at body 3. Case 2-1 is a typical case in which 
two opposing jets are colocated. In Cases 2-2 and 2-3, two opposing jets are 
not colocated. 

For a class of problems, such as Case 2-T, the modal equations become 

y. + = fan- * = 1 , • • ■ , n (2.3) 

where fa the i th modal gain and the scalar control input u is bounded as 

- 1 < u < 1 (2.4) 

The problem is to find the control input u(f) which minimizes the performance 
index 

J = /*' dt = t } 

subject to Eqs. (2.3) and (2.4), and given boundary conditions. The time- 
optimal bang-bang solution of this problem with the rest-to-rest maneuvering 
boundary conditions has, in most cases, (2n — 1) switches, and the solution is 
symmetric about the mid-maneuver time tj/2 [12]. That is, for a case with 
(2n — 1) switches, we have the symmetric switching pattern given as: 


— ^2n f 2n—j ! 


j = 1» • ' • ,n 


(2.5) 


C/j 
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where tj is the j th switching time and t 2n = tj- 

A bang-bang input with (2 n — 1) switches is expressed as 

u(t) = '£B } u 1 (t~t>) ( 2 - 6 ) 

:=o 

where B } is the magnitude of a unit step function u,(t) at t r This function 
can be characterized by its switch pattern as: 

B = { Bo, B\, B 2 , ■ • ■ , B 2n } 

T = { to, ti, < 2 , • • • , t 2n } 

where B represents a set of B : with Bo = B 2n = ±1 and B : = ±2 for j = 
1, • • • , 2n — 1; T represents a set of switching times (tj, • • • , < 2 n-i) and the initial 
and final times (to = 0 and tj — t 2n )- 

The rest-to-rest maneuvering constraint for the rigid-body mode (u^i = 
0) can be found as 

= 0 (2.7) 

z 1—0 

The i th structural mode solution for the control input of Eq. (2.6) for 

t > tj is 

<t> 

y.(0 = -^^BjCosuiiit-tj) 

u i J * 0 

d> 

= — — j[cosu>i(t — t n ) ^2 Bj cosui(tj — t n ) 

js =0 
2 n 

+ sinw,(t - t n ) Bj sin uii(tj - f n )] (2.8) 

j* 0 

and it can be shown that the following constraint equation for each mode 

2 n 

^Bj sinu^tj - t n ) = 0 
o 


(2.9) 
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is always satisfied for any bang-bang input which is symmetric about the mid- 
maneuver time t n . Consequently, we have the following flexible mode con- 
straints for no-residual structural vibration (i.e., y,(f) = 0 for t > tj ): 

3n 

£ Bj cos Ui(tj -t n ) = 0 (2.10) 

for each flexible mode. 


r* 

lUn, *1 


tj 

A u 2 or U 3 


Au- 


l N-l 



Figure 2.2: Pulse sequences. 

On the other hand, for Case 2-2 with the pulse sequences as illustrated 
in Fig. 2.2, the boundary conditions of the rest-to-rest maneuvering problem 
result in the following constraint: 

N-i n 

A i ~ ^ia 52 Aj =s 0 
i*o.a j«i4 


( 2 . 11 ) 
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where <j>\\ and o n are the modal input distribution coefficients associated with 
the rigid-body mode and the two control inputs u x and u 2 , and A ; and A' are 
defined as in Fig. 2.2. 

The positioning constraint for the rigid-body mode, with specified 
t/i(< > //), then becomes 

2y, (i/) - £ [2(A + Aj) (2.12) 

.7=0,2 

+ ^2 [2<; Aj + A*] = 0 

J = 1.3 

Also, rest-to-rest maneuvering requires that for each flexible mode, y,(t) = 0 
for t >t/\ i.e., we have 


N-l 

iV 


i 

+ 4>i1 = ® 

(2.13) 

J = 0,2 

]= 1.3 


N-l 

N 


s ij 

— 4>i7 ^2 S ij = 0 

(2.14) 

J = 0,2 

■ jisl,3 



where 

dj = cos(wifj) - cos(ui(tj + Aj)) 
Sij = sin(w,fj) - sin^t, + A,)) 


for i — 2, • • • , n. 

Remark: For cases in which the control inputs are one-sided, each 
control input for the time-optimal solution need not be an odd function about 
the mid-maneuver time. Thus, the problem with one-sided control inputs be- 
comes more difficult to solve them the standard problem with two-sided control 
inputs, and many theoretical issues (e.g., the uniqueness and structure of time- 
optimal solutions) need to be resolved. We now present a detailed solution of 


each case. 
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2.2.2 Case 2-1 with a Two-Sided Control Input 

In this case, as illustrated in Fig. 2.1, a two-sided control input is 
bounded as: 

-1 < u(t) < +1 

and the equations of motion for this case are 

m x xi + ki(xi — x 2 ) = u 

mjXj + &i(x2 — Xi) + k 2 (x 2 ~ x 3 ) = 0 

m 3 x 3 + k 2 (x 3 - x 2 ) = 0 

where Xj, X2 and 13 are the positions of body 1, body 2 and body 3, respec- 
tively, and the ’inal parameters are mj = m 2 = m 3 = ki = k 2 = 1 with 
appropriate unit., and time is in units of second. 

The boundary conditions for a rest-to-rest maneuver are given as 


ix(0) = x 2 (0) = x 3 (0) = 0 

(2.15) 

I 1 (*/) = *2(*/) = X3 (tf) = 1 

(2.16) 

ii(0) = x 2 (0) = x 3 (0) = 0 

(2.17) 

ii(tf) = x 2 (t } ) = x 3 (t f ) = 0 

(2.18) 

The modal equations are 

y x = 0.3333u 

(2.19) 

j/2 + <*>2 1/2 = 0.5u 

(2.20) 

i/3 + w 3 ya = — 0.1667u 

(2.21) 

u] = [-6- -Ak x k 2 c\/2 

(2.22) 


where 
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u\ = [-6 + \Jb 2 - ikiki c ]/2 
, —k\(m\ + m 2 )m 3 — A: 2 (m 2 + m-i)m\ 

b _ — — 

m 1 m 2 m 3 

C = (m 1 + m 2 + rn 3 )/(m l m 2 m 3 ) 

(2.23) 

(2.24) 

(2.25) 

and u >2 = 1 rad/sec and = \/3 rad/sec for the nominal system, 

sponding boundary conditions for the modal coordinates are 

The corre- 

yi( 0 ) = 5/2(0) = 5 / 3 ( 0 ) = 0 

(2.26) 

Vi(tf) = 1. 5/2 (</) = 5/a(*/) = 0 

(2.27) 

yi( 0 ) = 5 / 2 ( 0 ) = 2 / 3 ( 0 ) = o 

(2.28) 

= 5/2 (</) = y3(*/) = ° 

(2.29) 


The time-optimal control input with 5 switches is expressed as 

u(*) = u,(*) - 2u,(* - f : ) + 2u,(< - t 2 ) - 2u s (t - * 3 ) 

+ u,(t - *<) - u,(t - f 5 ) + u,(< - * 6 ) (2.30) 

with the symmetry conditions 

*6 = 2*3 
ts — 2*3 — ti 
ti = 2*3 — *2 

The time-optimal control problem is then formulated as the following 
constrained minimization problem: 

min J = tf = *$ = 2*3 (2.31) 

subject to the following constraints: 

6-*6 + D - 1 ) j + 1 [ 2 (*«-* 2) 2 ]=0 

j=i 


(2.32) 
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1 + cos(w 2 f 3 ) + 2 £(-l y cos(u; 2 (<3 - t>)) = 0 (2.33) 

;=i 

2 

1 + cos(u> 3 * 3 ) 4- 2 ^(-l)- 7 cos(u> 3 (f 3 - tj)) = 0 (2.34) 

;=i 

fi> h, <3 > 0 (2.35) 

A numerical technique of nonlinear constrained optimization problem discussed 
in Appendix A, was used to obtain a solution as: 

u(t) = u,(t) - 2u,{t - 0.944) + 2u a {t - 2.012) 

- 2u,(* - 3.255) + u,(t - 4.499) 

- u,(f- 5.567) + u,(<- 6.511) (2.36) 

The time responses of x 3 to this time-optimal control input axe shown in Fig. 2.3 
for four different values of A: = A:j = Ar 2 - We notice that the responses are quite 
sensitive to variations in the model parameters. Similar responses can also be 
observed for arbitrarily combined variations of A:, and m u but keeping the total 
mass constant (m j + rn 2 + m 3 = 3). For convenience, simulation results only 
for k = ki = k 2 variations are presented in this chapter. 

2.2.3 Case 2-2 with Two One-Sided Control Inputs 

For Case 2-2, as illustrated in Fig. 2.1, two one-sided control inputs 
axe bounded as 

0 < U! < +1 (2.37) 

-1 < u 2 < 0 (2.38) 

For this case, the time-optimal control inputs are assumed as: 

= u,(t) - u 3 (t - A 0 ) 
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+ u,(t - t 2 ) - u,(t - t 2 - £ 2 ) 

u 2 = -u,(t - ti) + U,{t - t l - Aj) 

(2.39) 


— U,(t — t 3 ) + U s (t — <3 — A 3 ) 

(2.40) 


The modal equations for this case are 


— 

jii = 0.3333ui + 0.3333u 2 

(2.41) 


y 2 +uj]y 2 = 0.5u! 

(2.42) 


$3 + uli /3 = 0.1667ui — 0.3333u 2 

(2-43) 

— 

and the rest-to-rest maneuver constraints are 


— 

A 0 — Ai + Aj — A 3 = 0 

(2.44) 


6 + D-imJ + 21^1 = 0 

J =0 

(2.45) 


51 [cos(u> 2 fj) - cos(u 2 {t, + Aj))] = 0 

J= 0,2 

(2.46) 


51 (sin(u> 2 *,) - sin(w 2 (t J + A,))] = 0 

:=Q,2 

2 53 [cos(watj) - cos(u^(fj + A,))] 

(2.47) 


J=1.3 



+ 5Z [co*(ws<i) ~ cosfwafc + Aj))] = 0 

J= 0,2 

(2.48) 


2 [sin(w 3 tj) - sin(u; 3 (tj + A,))] 

i*i.3 


— 

+ 51 («»(«•*<>) - sin(w 3 (t ; + A_,))] = 0 

J= 0,2 

(2.49) 

— 

The time-optimal solution can then be obtained by solving the constrained 
minimization problem: 


min J = tf = t 3 + A 3 

(2.50) 


subject to the constraint given by Eqs. (2.49). 
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A solution of this problem can be found as 

to = 0.0000, A 0 = 1.0429 

= 1.9581, Ai = 1.2401 

<2 = 3.1416, A 2 = 1.0429 (2.51) 

< 3 = 4.504, A 3 = 0.8457 
t } = 5.349 

The time responses of x 3 to the time-optimal control inputs are shown in 
Fig. 2.4 for four different values of k = k x = k 2 . Similar to Case 2-1, the 
responses are sensitive to variations in the model parameters. An interesting 
feature of this case is that the pulse sequences are of a “bang-off- bang” type, re- 
sulting in the control on-time of 4.172 sec, which is different from the maneuver 
time of 5.35 sec. 


2.2.4 Case 2-3 with Two One-Sided Control Inputs 

For Case 2-3, as illustrated in Fig. 2.1, two one-sided control inputs 
are bounded as 


0< <+l (2.52) 

-1 < u 3 <0 (2.53) 


Similar to Case 2-2, the time-optimal control inputs are assumed as: 


= «•(<) - - A 0 ) 


+ u,(t - t 2 ) - u,{t - t 3 - A 2 ) (2.54) 

“3 = -U,(<-*l) + «,(<-<! - Aj) 

- u s (t - h) + u,(t - 1 3-A3) (2.55) 


The modal equations for this case are 


yi = 0.3333(«i + U 3 ) 


(2.56) 


yi+uly? = 0 . 5 ( ui - u 3 ) ( 2 . 57 ) 

y3 + u>3j/ 3 = 0.1667(u, + u 3 ) (2.5S) 

and the rest-to-rest maneuver constraints are 

A 0 - A, + A 2 - A 3 = 0 (2.59) 

6 + E(— 1) J (A? + 2< j A j 2 ] = 0 (2.60) 

J=0 

3 

H(-l) J [sin(u/ 2 <j) - sin(u; 2 (t J + A,))] = 0 (2.61) 

;= o 

3 

^(-l) J [cos(w 2 t ; ) - cos(u> 2 (<j + Aj))] = 0 (2.62) 

:=o 

3 

^(-l) J [sin(u; 3 t ; ) - sin(u; 3 (< J + A i ))J = 0 (2.63) 

;=0 
3 

^(-lyieos^) - cos(u/ 3 (<_, + A,))] = 0 (2.64) 

:=o 

*i><2i*3><4,<5 > 0; <o = 0 


A : > 0 

The. time-optimai solution can then be obtained by solving the con- 
strained minimization problem: 

min J = t f = < 3 + A 3 (2.65) 

subject to the constraints given by Eqs. (2.64). 

A solution of this problem is 

<o = 0.0000, A 0 = 0.9510 
<i = 1.3631, Aj = 0.1658 

t 2 = 2.8329, A 2 = 0.1658 (2.66) 

t 3 = 3.4109, A 3 = 0.9510 
t f = 4.3619 

The time responses of z 3 to the time-optimal control inputs are shown 
in Fig. 2.5 for four different values of k = fcj = k 2 . Again, the responses are 



IS 


quite sensitive to variations in the model parameters. Similar to Case 2-2, an 
interesting feature of this case is that the pulse sequences are of a u bang-ofF- 
bang” type, resulting in the control on-time of 2.234 sec, which is different from 
the maneuver time of 4.362 sec. 

Compared to Case 2-1 and Case 2-2, this case has the fastest ma- 
neuver time as well as the smallest control on-time. Therefore, the actuator 
configuration of Case 2-3 can be considered to be “optimal” in the sense of 
minimizing both the maneuver time and the jet on-time. 

2.3 Robust Time-Optimal Control 

As shown in the preceding section, a standard, time-optimal control 
approach requires an accurate mathematical model and thus the resulting so- 
lution is often sensitive to plant modeling uncertainty. 

In this section, expanding on the approach introduced in Section 2.2, 
a parameter optimization problem is formulated with additional constraints for 
robustness with respect to the structural frequency uncertainty. The resulting 
robustified or desensitized , time-optimal soultion is a multi-switch bang-bang 
control, and is thus implementable for spacecraft equipped with on-off reaction 
jets [13]. 


2.3.1 Problem Formulation 


By taking the derivative of Eq. (2.8) with respect to u ly we get 


dyijt) 

du, 


cos Ui{t - ^-) Y^(tj - sin u t (t, - 


i-o 


(2.67) 
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for each flexible mode. Letting dy,(t)/duji = 0 for all t > tj, we have 

i = 2, ■ • • , n (2.68) 

J=0 

which are called the first-order robustness constraints for the case with a two- 
sided control input. 

Similarly, the robustness constraints for a case with two one-sided 
control inputs can be found as: 

N - 1 N 

~<t>, 1 E C *J + ^2 X] = 0 ( 2 - 69 ) 

j =0,2 j*l,3 

lV-1 S 

4>i\ s 'i ~ ^'2 E = ® (2.70) 

;=0,2 j=1.3 

where 


Cij = tj cos (Uitj) - (tj + A j) cos (u>i(tj + Aj)) 

5.-; = t : sin(u; t tj) - (tj + A ; ) sin + A,)) 

which are called the first-order robustness constraints for the case with two 
one-sided control inputs. 

2.3.2 Case 2-1 with a Two-Sided Control Input 

For Case 2-1, the time-optimal control is a five-switch bang-bang func- 
tion, but the resulting responses were shown to be very sensitive to variations 
in model parameters. A robustified, time-optimal solution of the same problem 
is now computed as follows. The robust time-optimal control input is assumed 
as: 

u (t) = u,(t) + 2 - *;)] + u,(t - ho) (2.71) 
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with the symmetry conditions 

t 6 = 2t$ — £4 

1 7 = 2t S — £3 

£§ — 2^5 — £3 (2.72) 

£9 = 2£s — t\ 
t\o = 2£s 

The constrained optimization problem with the first-order robustness 
constraint can be formulated as: 

min J = tj = £ 10 (2.73) 

subject to 

;=i 

4 

1 + cos(u> 2*5) + 2^[(-l) J cos(u > 2 (tj - *s))] = 0 

J=1 

4 

1 + cos(u> 3 f 5 ) 4- 2 £[(-l) J cos(u 3 (tj - * 5 ))] = 0 
j - 1 

4 

t s am(u 2 ti) + 2]£(-l) J (fj - f 5 )sin(w 2 (*j - <»)) = 0 
j=i 

4 

f 5 sin(w3t 5 ) + 2^(-l ) } (tj - t 5 )sin(u>3(t ; - f 5 )) = 0 

j=i 

*1>*2>*3»*4>*5 > <0 = 0 

A robust time-optimal solution with 9 switches can be found as: 

= 0.560, f 2 = 1.460 
t 3 = 2.690, U = 3.804 

*s = 5.091, t 6 = 6.377 (2.74) 

t 7 = 7.491, t 8 = 8.722 
t 9 = 9.622, tio = 10.18 

The time responses of x 3 to this robust time-optimal control input are shown 
in Fig. 2.6 for four different values of k = fcj = fc 2 . We notice that the re- 
sulting responses are less sensitive to parameter variations, as compared to the 
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responses to the ideal, time-optimal control input, as shown in Fig. 2.3. The 
second flexible mode is significantly excited during maneuvers, however. Per- 
formance robustness has been increased at the expense of the increased maneu- 
vering time of 10.18 sec, as compared to the ideal minimum-time of 6.511 sec. 
It is, however, emphasized that simply prolonging the maneuver time does not 
help to reduce residual structural vibrations caused by modeling uncertainty; 
a proper pulse sequence is necessary. 

2.3.3 Case 2-2 with Two One-Sided Control Inputs 

For Case 2-2, we can represent the control inputs as follows: 

uj = u 3 (t) - u s (t - A 0 ) + u,(t - t 2 ) 

- u s (t -t 2 - A 2 ) + u,(t - t<) - u s (t - t 4 - A 4 ) 

“2 = ““»(* ~ h) + - *1 - A!) - u s (t - t 3 ) 

+ u,(t -t 3 - A 3 ) - u,(t - *s) + u,(t -t s - As) 

where we have 11 unknowns to be determined, and tj and A ; are defined as in 
Fig. 2.2. 

The robust time-optimal solution can then be obtained by solving the 
constrained parameter optimization problem 

min J = tj — ti 4- A* (2.75) 

subject to 

Ao — Ai + Aj — A 3 + A 4 — Aj = 0 

6 + £(-im? + 2f,Aj] = 0 

j=0 

£ [cos(u> 2 tj) - cos (u 2 {tj + Aj))] = 0 

j= 0,2.4 
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Y [ain(wjt,-) - sin (u 2 (tj + Ay))] = 0 

js 0,2,4 

2 Y [cos(w 3 t ; ) - cos(w 3 (< ; + A ; ))] 

J = 1 .3 

+ 5Z ( cos ( w 3*j) - cos(w 3 (tj + Aj))] = 0 

2=0, 2, 4 

2 ^ [sin(w 3 ^) - sin(u> 3 (tj + A,))] 

;=i.3 

+ Y l«n(o*tj) - sin(u> 3 (tj + A,))] = 0 

;=0,2,4 

Y [tj CO s(u> 2 tj) - (tj + Aj) cos (v 2 (tj + Aj))] = 0 

i= 0,2,4 

r [*j sin(u> 2 fj) — ( t 3 + Aj) sin(u> 2 (fj + A 7 ))] = 0 

j=0,2,4 

2 H (tjCOs(w 3 tj) - (tj + Aj)cos(w 3 (fj + Aj))] 

J=l,3 

+ H [<; cos(u> 3 <j) - (tj + Aj) cos(w 3 (tj + Aj))] = 0 

j=0,2,4 

2 Y l*i sin(w 3 tj) - (tj + Aj) sin(u; 3 (t J + Aj))] 

i=i,3 

+ H (*j sin(w 3 tj) - (tj + Aj) sin(w 3 (fj + A ; ))] = 0 

j=0,2,4 ■ 

Aj > 0; j -0,1, 2, 3, 4, 5 

ti, t 3 , t 3 , t^, t$ ^ 0 ) to — 0 

A solution to this problem is: 

to = 0.0000, A 0 = 0.5814 
ti = 1.9948, Ax = 0.5657 
t 2 = 2.8533, A 2 = 1.2202 

t 3 = 4.3234, A 3 = 1.4425 (2.76) 

U = 6.3455, A< = 0.5814 
t 5 = 6.9850, A s = 0.3747 
tj = 7.36 

The time responses of x 3 to the robust time-optimal control inputs are 
shown in Fig. 2.7 for four different values of k = k x = k 2 . Similar to Case 2-1 of 
the preceding section, the robustness has been increased at the expense of the 
increased maneuvering time of 7.36 sec, as compared to the ideal minimum- 
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time of 5.35 sec. The jet on-time is 4.765 sec. The second flexible mode is less 
excited, as compared to Case 2 - 1 . 

2.3.4 Case 2-3 with Two One-Sided Control Inputs 

Assuming that each control input has two pulses as in Case 2 - 2 , we 
can represent the control inputs as: 

u x = u,(t) - u a (t - A 0 ) 4- u,(t - t 2 ) - u s (t - t 2 - A 2 ) 

+ u,(t - t 4 ) - u,(t - t 4 - A„) 
u 3 = -u 4 (f - ti) + u a (t - ti - Aj) - u,{t - t 3 ) 

+ — <3 — A 3 ) — U s (t — <5) + U s (t — <5 — A5) 

where we have 11 unknowns to be determined, and t : and A ; are defined as in 
Fig. 2.2. 

The robust time-optimal solution can then be obtained by solving the 
constrained parameter optimization problem 

min J = t 5 + A s (2-77) 

subject to 

Ao — Aj + A 3 — A 3 + A 4 — As — 0 

« + B-im? + 2‘A ! ] = o 

;=o 

5 

^(-lyjsin^tj) - sin (w 2 (t, + A,))] = 0 

j=0 

5 

52(-l) J [cos(u;2t J ) - cos(u 2 (tj + A,-))] = 0 
1=0 

O -1 ) - ' sin(w 3 tj) - sin (ua(t> + A,))] = 0 
;=o 
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0-1)' cos ( w 3 ^) - cos(u > 3 {t : + A,))] = 0 

J=0 

5 

cos U> 2 1 3 - (tj + Aj) COSU^f, + A ; )] = 0 
;=o 

5 

- (tj + AJsinu ) 2 (tj + Aj)j = 0 

i=o 

5 

cos Ufilj - (tj + A : )cosu %(tj + A ; )] = 0 

>* 0 
5 

““ (tj + Aj)sinw 3 (^ + A ; )] = 0 

j=o 

A, >0; j = 0, 1,2, 3,4, 5 

^3> ^4i ^5 > 0; to = 0 

A solution to this problem is 

to = 0.000, A 0 = 0.2189 
h = 3.774, Ai = 0.2609 
t 2 = 2.030, A 2 = 0.3594 

t 3 = 5.778, A 3 = 0.3594 (2.78) 

U = 4.152, A< = 0.2609 
t s = 7.968, A 5 = 0.2189 
tj = 8.187 

The time responses of z 3 to the robust time-optimal control inputs are 
shown in Fig. 2.8 for four different values oi k = k x = k 2 . The robustness has 
been increased at the expense of the increased maneuvering time of 8.187 sec, as 
compared to the ideal minimum-time of 4.362 sec. However, the control on-time 
is only 1.678 seconds! Because of the properly coordinated pulse sequences, the 
flexible modes are not significantly excited during maneuvers and the residual 
responses after the maneuvers are well desensitized. 

Now consider the effects of the accuracy of the switching time of firing 
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of the pulses. The solution of Case 2-3 is now truncated to 2 decimal places : 

to = 0.00, A 0 = 0.21 
t, = 3.77, A! = 0.26 
t 2 = 2.03, A 2 = 0.35 

t 3 = 5.77, A 3 = 0.35 (2.79) 

U = 4.15, A< = 0.26 
< 5 = 7.96, A s = 0.21 

tj = 8.18 

The time responses of x 3 to the robust time-optimal control inputs are shown 
in Fig. 2.9 for four different values of k = k x = k 2 . The residual vibrations have 
not increased and the response is robust to the switching time. Now consider 
truncating the solution to one decimal place, which may be too unrealistic in 
practice, as follows: 

to = 0.0, A 0 = 0.2 
fi = 3.7, Ai = 0.2 
t 2 — 2.0, A] — 0.3 

t 3 = 5.7, A 3 = 0.3 (2.80) 

U = 4.1, A« = 0.2 
t 5 = 7.9, A s = 0.2 

tj = 8.1 

The time responses of x 3 to the robust time-optimal control inputs are shown in 
Fig. 2.10 for four different values of k = k x = k 2 . Again the residual vibrations 
have not increased but the steady-state solution has an offset. 


2.4 Summary 

In this chapter, we have demonstrated that the proposed robustifica- 
tion or desensitization approach does generate robust time-optimal open-loop 
control inputs for an uncertain dynamical system. Furthermore, on the con- 
trary to a common notion, the results of this chapter indicate that properly 
coordinated, on-off pulse sequences can achieve a fast maneuvering time with 
a minimum of structural vibrations during and/or after a maneuver, even in 
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the face of plant modeling uncertainty. The time-optimal responses have been 
desensitized at the expense of the increased maneuvering time. It is again 
emphasized that simply prolonging the maneuver time does not help to re- 
duce residual structural vibrations caused by modeling uncertainty; a proper 
coordination of pulse sequences is necessary, as demonstrated in this chapter. 


Table 1: Summary of the results 


Time-Optimal Control 



Jet On-Time 

Case 2-1 

6.511 

6.511 

Case 2-2 

5.340 

4.172 

Case 2-3 

4.362 

2.234 

Robust Time-Optimal Control 


J m = tf (sec) 

Jet On-Time 

Case 2- 1 

10.18 

10.18 

Case 2-2 

7.360 

4.765 

Case 2-3 

8.187 

1.678 


The results of this chapter are summarized in Table 1. As can be 
noticed in this table, it is natural to select the actuator configuration of Case 
2-3, since this case provides the “best” overall performance in the sense of min- 
imizing the maneuvering time, fuel consumption (jet on-time), and structural 
mode excitation. For Case 2-1, the maneuvering time and the jet on-time are 
the same, which is clearly undesirable from the viewpoint of fuel consumption. 
To avoid such undesirable continuous jet firings during a maneuver, a robust 
fuel- and time-optimal control problem will be treated in the following chapter. 

2.5 Conclusions 

A time-optimal open-loop control problem of flexible spacecraft in 
the face of modeling uncertainty has been investigated. The primary study 
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objective was to explore the feasibility of computing open-loop, on-off pulse 
control logic for uncertain flexible spacecraft. The results indicate that the pro- 
posed approach significantly reduces the residual structural vibrations caused 
by modeling uncertainty. The results also indicate the importance of a proper 
jet placement for practical trade-offs among the maneuvering time, fuel con- 
sumption, and performance robustness. 
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Time (sec) 

Figure 2.4: Responses of time-optimal control Case 2-2. 
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x3 



Time (sec) 

Figure 2.6: Responses of robust time-optimal control Case 2-1. 

















Chapter 3 


Robust Fuel- and Time-Optimal Control of Uncertain 

Flexible Spacecraft 


3.1 Introduction 

The fuel- and time-optimal control problem of flexible spacecraft is 
in general a difficult problem to solve, even for a case without plant modeling 
uncertainty (e.g., see [14,15,16]). A standard, optimal control approach to such 
a problem requires an accurate plant model, and thus the resulting solution 
becomes sensitive to variations in model parameters. For this reason, a pulse- 
modulated, classical feedback controller is often employed for real spacecraft 
equipped with reaction jets [13,17]. 

In this chapter, however, we attempt to develop robustified, open- 
loop, fuel- and time-optimal control inputs for uncertain flexible spacecraft, 
which are often required to maneuver as quickly as possible without significant 
structural vibrations during and/or after a maneuver. In particular, we inves- 
tigate the rest-to-rest maneuvering problem of a flexible spacecraft equipped 
with on-off reaction jets in the presence of uncertainty in model parameters. 
Expanding on the recent results of [11,12], we formulate a constrained optimiza- 
tion problem, where the objective function to be minimized is a weighted sum 
of the consumed fuel and the maneuvering time, with additional constraints 
for robustness with respect to plant modeling uncertainty. 

We assume that the mass of the consumed fuel is small compared 
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with the mass of the spacecraft, and that the consumed fuel is proportional to 
the jet on-time. We further assume that the structural flexibility and mass dis- 
tribution of the vehicle are more uncertain than the total mass of the system. 
Consequently, we focus on the robust control problem of flexible spacecraft in 
the face of modal frequency uncertainty. Furthermore, many theoretical and 
practical implementation issues inherent to a constrained parameter optimiza- 
tion problem are not elaborated in this chapter. The ultimate research goal of 
the future is to develop nonlinear feedback control logic for achieving the robust 
time-optimal performance similar to that presented in this chapter. As a first 
step toward such a research goal we emphasize the formulation of a robust fuel- 
and time-optimal control problem of flexible spacecraft and we investigate its 
solution in terms of switching patterns. 

Other robustified, open-loop approaches attempt to find a smooth 
continuous forcing function (e.g., a versine function) that begins and ends with 
zero slope. The basic idea behind such approaches is that a smooth control 
input without sharp transitions is less likely to excite structural modes during 
maneuvers. The results of this chapter further confirm that most open-loop 
approaches, which utilize a smooth continuous control input so that structural 
modes are less likely to be excited, do not fully utilize the available control 
energy in performing a time-optimal maneuver. 

The remainder of this chapter is organized as follows. Section 3.2 de- 
scribes the fuel- and time-optimal control problem of flexible spacecraft without 
modeling uncertainty. The control problem is transformed into a parameter op- 
timization problem in which the objective function is a sum of the final time 
and the jet on-time. A simple dynamical system with a rigid-mode and one/two 
flexible mode, shown in Fig. 3. 1-3.3, is used to illustrate the concept. Three 
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cases ( Case 3-1, Case 3-2 and Case 2-1) are explored (see Figs. 3. 1-3.3). Sec- 
tion 3.3 describes the same problem that of Section 3.2 but with consideration 
of the presence of modeling uncertainty. Robustified fuel- and time-optimal 
solutions are then compared with an ideal solution obtained in Section 3.2. 

3.2 Fuel- and Time-Optimal Control 

3.2.1 Problem Formulation 

Consider a flexible spacecraft described by 

Mx -f Kx — Gu (3.1) 

where x is a generalized displacement vector, M a mass matrix, K a stiffness 
matrix, G the control input distribution matrix, and u the control input vector. 

Equation (3.1) is transformed into the modal equation: 

y + Q 7 y = (3.2) 

where y is the modal coordinate vector, f I 7 = diag( u>f), the i th modal fre- 
quency (u>! = 0 for the rigid body mode), and $ the modal input distribution 
matrix. 

In this chapter, we consider a simple model of a flexible spacecraft 
with a rigid-body mode and one dominant flexible mode, as shown in Fig. 3.1. 
We explore two cases for this system with two control inputs, u + and u_ , which 
are bounded as 

0< u+ <+l (3.3) 

-1 < u_ <0 (3.4) 

For Case 3-1, both u + and u_ are acting on body 1, resulting in a typical case 
with a two-sided control input u\ with |ui| < 1. For Case 3-2, u + is placed 
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on body 1 and u_ is placed on body 2, resulting in a case with two one-sided 
control inputs: u x = u + and u 2 = u_. For Case 2-1 both u+ and u_ are acting 
on body 1, of the three mass-spring system. 

The boundary conditions for the rest-to-rest maneuver are expressed 
in modal coordinates as: 

yi(0) = 0, yi(t/) = 1 

y i(0) = 0, t/i(</) = 0 . . 

y,(0) = o, yi {t } ) = o 1 oj 

y.(o) = o, t )i{t } ) = o 

for i = 2, n 

For a fuel- and time-optimal control problem, the objective function 
to be minimized is, in general, a weighted sum of the maneuvering time and 
the consumed fuel. However, we consider here the following objective function 
which is simply a sum of the maneuvering time and the product of weight a 
and total jet on-time: 

J = f 7 {1 + a(|u+| + |u_|)}d< (3.6) 

Jo 

The problem is then to find the control inputs which minimize the performance 
index J subject to Eqs.( 3.2)-( 3.4). In Section 3.3, we investigate the same 
problem as above but with additional robustness constraints with respect to 
plant modeling uncertainty. 


3.2.2 Rigid Body Fuel- and Time Optimal Control 

For a rigidized model of the nominal model discussed the equation of 

motion is 


(mi -I- m 2 )z = u 


(3.7) 
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The rest-to-rest manuever, fuel- and time-optimal solution for the objective 
function given by Eq.(3.6) can be solved using control input as a combination 
of unit-step function u s as 

u = u,(t) - u s (t - A) - u,(t - t x ) + u,(t -ti- A) (3.S) 

The solution is A = j(m x + m 2 )/( 2a + 1) and f, = 2A(a + 1) = 3.27 sec 
for a = 1. When this solution is applied to the flexible model there is residual 
vibration present and it remains even by prolonging maneuvering time for large 
weight a. 

3.2.3 Parameter Optimization Problem 

The control inputs u+ and tt_ can be expressed as a combination of 
unit-step function u, as 

JV-l 

= E M* “ l i) ~ “•(* -*j- A>)] (3.9) 

>* 0,2 

jV 

u - = 5Z [-«,(<- t>) + «,(<- <> - A>)] (3.10) 

>= 1.3 

which are in the form of one-sided pulse sequences as shown in Fig. 3.4. The 
j pulse starts at tj and ends at ( t } + A^). A j is the j til pulse duration. Due 
to the symmetric nature of the rest-to-rest maneuvering problem, we assume 
that u + and u_ have the same number of pulses, (iV+ l)/2, where N is defined 
as in Fig. 3.4. 

Substituting Eq. (3.10) into Eq. (3.2) and incorporating the boundary 

conditions of Eq. (3.5), we get the time response of the rigid-body mode as 

, N-l 

Vi(< > </) = ~\<hi E - 2tjAj - A?) 

1 j=t0,2 

-*12 E (2tA, - 2tjAj - AJ) ] 

>* 1.3 


(3.11) 
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where 0, ; is the (i,j) element of the modal input distribution matrix $. For a 
given boundary condition, y\{t > tj) = 1, the following constraint must hold 

n-\ N 

4> 11 2 E a j = 0 (3.12) 

2=0,2 > = 1.3 

The rest-to-rest constraint for the rigid-body mode then becomes 
2 — 4>\\ ^2 (2 tj&j — A*] + 22 [2 tj&j ~ ^;] = 0 (3.13) 

2=0,2 2=1.3 

Substituting Eq. (3.10) into Eq. (3.2) and solving for the time response 
of the i th flexible mode, we get 

l N-l N 

Vi(t) = — COS (w.Oh^i 22 C 'j + 4,2 J2 c 2 

u> ' 2=0.2 >=1,3 

j N-l N 

"b j sin(u2|i)[0,'i ) ) Sij <f >{2 ) ' -s,>] (3.14) 

U ' 2=0,2 >=1,3 

where 


Cij = cos (u,tj) - cos(u Ji(tj + Aj)) 

Sij = sin(u>,fj) - sin(u/,(tj 4- A,)) 

for t >tf. 

Also, rest-to-rest maneuvering requires yi(t) = 0 for t > tj\ i.e., we 
have the following flexible mode constraints for no residual structural vibration: 


jV-1 N 

-<t>i i co + ^*2 Y1 c o = ° 

(3.15) 

;=0,2 ; = 1,3 

JV-1 AT 


<f>n ^2 $ij - <i>i7 Yi s o = 0 

j— 0,2 2*1,3 

(3.16) 


for each flexible mode. 
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Since Eqs. (3.12), (3.13), and (2.14) are the necessary conditions for 
optimal solution, the fuel- and time-optimal control problem can be formulated 
as a constrained parameter optimization problem as follows: 

;V 

min J = tj + a £ A, ( 3 . 17 ) 

j— 0,1,2 

subject to the constraints given by Eqs. (3.12), (3.13) and (3.16). 

3.2.4 Case 3-1 

Consider a case where two control inputs u+ and u_ both are acting 
on body 1. The control inputs are then assumed as: 

u+ = u s (t) - u t (t - A) (3.18) 

u_ = -u,(t - h) + u,(t — f, — A) (3.19) 

where the pulse sequence is defined as shown in Fig. 3.2, and the maneuver time 
tf = i i + A. The nominal parameters are: TOj = m? — k = 1 with appropriate 
units. The corresponding matrices S7 2 and $ in Eq. (3.2) are 


■ 0 0 ' 

^ 1 

■ 1 1 ■ 

0 2 

* = 2 

1 1 j 


The constraint equations can be obtained as : 

t/A - A 2 - 2 = 0 (3.20) 

cos(u>f;/2) - cos(u(tf/2 - A)) = 0 (3.21) 

where u> = >/2. 

The fuel- and time-optimal solution is obtained by solving the con- 
strained minimization problem: 


min J — tj + 2A 
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subject to the preceding constraints. The solution of this problem is obtained 
as 


A = 0.45016 

tj = 4.44288 

t f = 4.89304 (3.22) 

The time responses of x 2 (the position of body 2) to the fuel- and 
time-optimal control inputs are shown in Fig. 3.5 for four different values of 
spring stiffness k. For the nominal case with k = 1 , the jet on-time is now 
significantly reduced by 79% while the maneuvering time is increased by 18%, 
comparing to the time-optimal case. (The time-optimal solution for Case 2- 
1 of chapter 2 has both maneuvering time and control on-time of 4.128 sec.) 
As expected, however, significant residual structural vibrations can be seen in 
Fig. 3.5 due to variations in the uncertain parameter k. 

3.2.5 Case 3-lb 

Consider a special case of identical to Case 3-1 but with a longer 
maneuver time by increasing the weight a on the fuel in the cost function. 

min J = tf + 2a A 

subject to the constraints given by Eqs(3.21). The solution of this problem for 
a = 10 is obtained as 


A = 0.225 
t\ — 8.88 
tj = 


9.11 


(3.23) 
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The time responses of x, (the position of body 2) to the fuel- and time- 
optima! control inputs are shown in Fig. 3.6 for four different values of spring 
stiffness k. The residual vibrations are still significant and are not reduced 
by increasing the maneuver time. Robustness to parameter variations is not 
achieved by increasing maneuver time. 


3.2.6 Case 3-2 

Consider Case 3-2 shown in Fig. 3.2 with two control inputs u, and 

acting on body 1 and body 2, respectively. The control inputs are assumed 
as 


“I = u s (t) - U 3 (t - A) 

u 2 = -u s (t - <j) + u s (t - tj - A) 

and the maneuver time *,«<, + A. The matrices & and $ are 
The constraint equations can be obtained as : 


Q 2 = 


0 0 
0 2 


(3.24) 

(3.25) 


tj - 2/A - A = 0 

sin(u>t//2) + sin(w(A - tj/2)) = 0 


(3.26) 

(3.27) 


The fuel, and time-optimal solution is obtained by solving the constrained 
nunirruzation problem: 


min J = t/ + 2A 


subject to the preceding constraints. 
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The solution is obtained as 

A = 0.9003 
= 2.2215 
t } = 3.1218 

The time responses of x 2 to the fuel- and time-optimal control inputs 
are shown in Fig. 3.7 for four different values of spring stiffness k. Comparing 
to the time optimal case, the jet on-time is now significantly reduced by 56% 
while the maneuvering time is also reduced by 24%. This result indicates the 
importance of proper jet placement in the fuel- and time-optimal problem. As 
expected, significant residual structural vibrations can be seen in Fig. 3.7 due 
to variations in the uncertain parameter k. 


3.2.7 Case 2-1 

For Case 2-1, as illustrated in Fig. 3.3, two one-sided control inputs 
are bounded as 


0 < u + < +1 ■ 

(3.28) 

O 

VI 

1 

VI 

H 

1 

(3.29) 


For this case, the time-optimal control inputs are assumed as: 

U! = u,(£) - u s (t - A 0 ) 

■f u*(£ — £j) — Uj(t — £2 — A 2 ) (3.30) 

Uj = — u,(£ — £1) + u*(£ — £1 — Aj) 

— u,(t — £3) + u,(£ — £3 — A3) 


(3.31) 
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The modal equations for this case are 


yx = 0.3333(u+ -f u_) 

(3.32) 

— 

h + ^ly7 = 0.5 (u+ + u_) 

(3.33) 


$3 + 1^3 2/3 — 0.1667(u+ + u_) 

(3.34) 

*—* 

where = 1 and = \/3. The rest-to-rest maneuver constraints are 


— 

A 0 — Ai + Aj — A 3 = 0 

3 

(3.35) 


6 + D-l) J [A?-2t J A i ] = 0 
i® 0 

(3.36) 


3 

EC-D^W^t,) - cos (uj 2 (tj + Aj))] = 0 

(3.37) 

— 

i®*o 
. 3 

2^(-l) J [sin(cj 2 *;) - sin(w 2 (< ; + A,-))] = 0 
;=o 

(3.38) 


3 

^(-l) J [ cos (w 3 t ; ) - cos(u> 3 (* ; + Aj))] = 0 
0 

(3.39) 


3 

^(-l)' ? [sin(w 3 ^) - sin(w 3 (ij + Aj))] = 0 
;=0 

(3.40) 

— 

The time-optimal solution can then be obtained by solving the constrained 


minimization problem: 

*5 

sJ 

min J = if + ^2 Ai 

(3.41) 



i=0 

subject to the constraint given by Eqs. (3.40). 


A solution of this problem can be found as 

to = 0.0000, A 0 = 0.400 

= 4.536, Ai = 0.135 

t 2 = 2.582, A 2 = 0.135 (3.42) 

t 3 = 6.854, A 3 = 0.400 
if = 7.251 

The time responses of x 3 to the fuel and time-optimal control inputs are shown 
in Fig. 3.8 for four different values of k = k x = k 2 . Now consider the time 
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response of the rigid body fuel- and time-optimal control solutions with final 
time tj = 7.35 seconds. The residual vibrations are much larger and the flexible 
mode is significantly excited as shown in Fig. 3.9.(Rigidized Case 2-1) 

3.2.8 Case 2-lb 

Consider a special case of identical to Case 2-1 but with a longer 
maneuver time by increasing the weight a on the fuel in the cost function. 


mm 


J = t f + a A t 


(3.43) 


1=0 


subject to the constraint given by Eqs. (3.40). 

A solution of this problem for a = 14 can be found as 

to = 0.0000, A 0 = 0.153 
U m 11.762, A x = 0.092 

t 2 = 2.152, A 2 = 0.092 (3.44) 

t 3 = 13.854, A 3 = 0.153 
t f = 14 

The time responses of x 3 to the fuel and time-optimal control inputs are shown 
in Fig. 3.10 for four different values of k = ki = k 2 . The residual vibrations 
are still significant and are not reduced by increasing the maneuver time. Ro- 
bustness to parameter variations is not achieved by increasing maneuver time. 


3.3 Robust Fuel- and Time-Optimal Control 

3.3.1 Problem Formulation 

As illustrated in Figs. 3.5 and 3.7, a standard, optimal control ap- 
proach requires an accurate plant model and thus the resulting solution is not 
robust to plant modeling uncertainty. For this reason, an open-loop optimal 
control approach is seldom used in practice. In this section, we attempt to 
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develop robust, open-loop, fuel- and time-optimal control inputs for flexible 
spacecraft. Expanding on the results of the previous section, we formulate 
a parameter optimization problem, where the objective function to be min- 
imized is a weighted sum of the consumed fuel and the maneuvering time, 
with additional constraints for robustness with respect to structural frequency 
uncertainty. 

To attenuate residual vibrations of the flexible modes, the energy of 
the residual vibrations should be minimized, where the energy of the residual 
vibration is proportional to the square of the amplitude. Thus, Eq. (3.14) can 
be written as 


Vi(t > tj) = Asin(u;;f) + f? cos(u>,f) (3.45) 

where A, B are functions of w t - and 4>{j and 


C 2 = A 2 + B 2 (3.46) 

dC 2 dC 2 du k 

~J~ = 1 T~ (3-47) 

dp, du k dp, 

where p, is the i th uncertain parameter. For these derivatives to be zero ( t > t f ), 
we get 


where 


jV-l n 

^2 9ij + 4>i2 £ 9ij = 0 

>=0,2 >=1,3 

N - 1 N 

<t>n £ hi] - 4>ii £ h,j = 0 

J=0,2 >=1,3 


(3.48) 

(3.49) 


9ij = ^ cos(u,tj) - (tj + Aj) 


cos(u>i(tj + Aj)) 


(3.50) 


hij = tj sin(w,tj) - (tj + Aj) 
sin (u,(t 3 + Aj)) 


(3.51) 
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Eqs. (3.48)-(3.49) are called the first-order robustness constraints. Similarly, 
the reorder robustness constraints can be expressed as 


{t,) m cos(u/.-tj) 


-(tj + A ; ) m cos (Wi(i> + A,)) 

(3.52) 

(ij) m sinMj) 


— (tj + Aj) sin(u>,(0 -f- Aj)) 

(3.53) 


for m = 1 , 2 , • • • , r. 

These robustness constraints are incorporated in the constrainted pa- 
rameter optimization problem formulation. Consequently, the number of pulses 
for each control input is changed to match the increased number of constraints. 
As an example, we consider the first-order robustness constraint, incorporated 
with the rest-to-rest maneuver constraints, to determine robust fuel- and time- 
optimal pulse sequences. ■ 


3.3.2 Case 3-1 

Assuming that each input has two pulses, we can represent the control 
inputs as follows: 


MO - u *(t - A 0 ) + U,(t - t 2 ) 


A 2 ) 

(3.54) 

-u s (t - ti) -1- u,(t - ti - AO - u,(f - t 3 ) 


+u,(t — t 3 — A 3 ) 

(3.55) 


where there are 7 unknowns to be determined, and tj and A : are defined as 

shown in Fig. 3.3. 
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The robust fuel- and time-optimal solution can then be obtained by 
solving the constrained parameter optimization problem: 

min J = £3 + A3 + (Aq + Ai 4- Aj + A3) ( 3 . 56 ) 


subject to 


Ao — Ax + A? — A3 — 0 

£(-l)’[2i,A, + Aj] + 4 = 0 

3=0 

3 

— lyicoswfj — cos u(t : + Aj)] = 0 
3 = 0 
3 

— sinu)(^ + A;)] = 0 

3=0 

3 

cos utj - (tj + Aj)cosw(fj + Aj)] = 0 

j=o 

3 

sinw^ - (tj + Aj)sinw(*j + A ; )] = 0 

3=0 

A, > 0; j = 0,1,2, 3 

tu *2> h > 0 


A solution to this problem is obtained as: 
to = 0.0000, A 0 = 0.2379 

U = 4.1575, Aj = 0.2489 . . 

t 2 = 2.3732, A 2 = 0.2489 } 

t 3 = 6.5418, A 3 = 0.2379 

The time responses of x 2 to this robust fuel- and time-optimal control 
inputs are shown in Fig. 3.11 for four different values of k. It can be seen that 
the robustness has been significantly increased at the expense of the increased 
maneuvering time of 6.78 sec. The jet on-time is 0.97 sec. Fig. 3.11 illustrates 
the reduction in residual vibration. 
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3.3.3 Case 3-2 

Now consider Case 3-2 with two control inputs expressed as: 

ui = u s (t) - u s (t - A 0 ) + t2,(* - t 2 ) 

— u s (t — t 2 — A 2 ) 

U 2 = -u,(< -<l)+ u,{t - ti - A!) - U s (t - t 3 ) 

+ U,(t — tj — A3) 

where there are 7 unknowns to be determined. 

The problem is to minimize the following objective function 

min J = + A3 + (Ao + Ai + A2 + A3) 

subject to 

Ao ~ Ai + A2 ~ A3 = 0 

i;(-l) ) [2l J A J +Aj]+4 = 0 

1=0 

3 

^[coswtj — cos ui(tj + Aj)] = 0 

j=o 

3 

^[sinu>^ — sinu>(tj -f A ; )] = 0 

j=o 

3 

y>, cos uitj — ( tj + Aj) costj(/j + Aj)] = 0 
2=0 
3 

y.ft, sinurtj — (tj + A ; )sinu/(tj + Aj)] = 0 
2=0 

Aj > 0; j = 0,1,2, 3 


(3.58) 

(3.59) 

(3.60) 


*1) *2) *3 > 0 
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A solution to this problem is obtained as: 

t 0 = 0.0000, A 0 = 0.3762 

= 2.5918, Ai = 0.4238 ngn 

£ 2 = 1.9151, A 2 = 0.4238 1 ' ’ 

<3 = 4.5545, A 3 = 0.3762 

The time responses of to this robust fuel- and time-optimal control 
inputs are shown in Fig. 3.12 for four different values of k. 

From Fig. 3.12, we notice that the residual vibration has been reduced 
considerably. 

3.3.4 Case 2-1 

For Case 2-1, we can represent the control inputs as follows: 

u+ = Uj(£) — Uj(t A 0 ) + u,(£ £2) 

— u,(t - t 2 - A 2 ) + u,{t - f<) - u s (t - U - A 4 ) 

u_ = —u,(t — t\) + u,(£ — — Aj) — u,(£ — £3) 

+ u,(£ — £3 — A3) — Uj(£ — £5) + u s (t — £5 — A5) 

where we have 11 unknowns to be determined, and t 3 and A ; are defined as in 
Fig. 3 . 3 . 

The robust time-optimal solution can then be obtained by solving the 

constrained parameter optimization problem 

s 

min J = tj 4- ^ A,- (3.62) 

.=0 


Ao — A] + A 2 — A 3 + A 4 — A5 — 0 

6 + y;(-i)'[a ) I -2< i a i ] = 0 

1=0 


subject to 


53 


£(-l) J [cos(u^) - cos(w 2 (<; + A,-))] = 0 
;=o 

^(-l) ; [sin(w 2 *j) - sin(w 2 (ij + A,-))] = 0 

;=o 

^(-l) J [cos(w 3 *j) - cos (w 3 (*j + Aj))] = 0 

;=0 

— l)’ ? fsin(w 3 ^j) - sin(u> 3 (*j + Aj))] = 0 
2 =o 

5 

^(-1)-% cos(w 2 /j ) - (tj + A ; )cos(u > 2 {tj + Aj))] = 0 . 

;=0 

£(-l) J [*., sin(u> 2 *,) - {tj + A J -)sin(wa(i > + A,))] = 0 

j=° 

£(-l) J ’[fjCOs(u> 3 *j) - {t } + Aj)cos (w 3 (fj + Aj))] = 0 
;=0 

^(-l) J [/j sin(u> 3 fj) - (ij + Aj)sin(w 3 (*j + Aj))] = 0 
; = 0 

Aj > 0; j -0,1, 2,3, 4, 5 

*1, *2i *3» *4 > *5 > 0 ! *o = 0 

A solution to this problem is: 

*o = 0.0000, A 0 = 0.1209 
*,=6.6158, A, = 0.1411 
* 2 = 2.3022, A 2 = 0.2018 

* 3 = 9.0502, A 3 = 0.2018 (3.63) 

t< = 4.7973, A„ = 0.1411 
* 5 = 1L433, A s = 0.1209 
t f = 11.554 

The time responses of x 3 to the robust fuel and time-optimal control 
inputs are shown in Fig. 3.13 for four different values of k = k\ = k 2 . The jet 
on-time is 0.927 sec. Now consider the time response of the rigid body fuel- and 
time-optimal control solutions with final time */ = 11.2 seconds. The residual 
vibrations are much larger and the flexible mode is easily excited as shown in 
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Fig. 3.14. Thus the flexible modes cannot be neglected in the model for they 
can be easily excited. 


Now consider the effects of uncertainty in the switching time of firing 
of the pulses. The solution of Case 2-1 is now truncated to 2 decimal places : 


to = 0.00, 

fi = 6.61, 

f 2 = 2.30, 
f 3 = 9.05, 
U = 4.79, 
f 5 = 11.4, 
tj = 11.5 

The time responses of x 3 to the robust i 


A 0 = 0.12 
A, = 0.14 
A 2 = 0.20 

A 3 = 0.20 (3.64) 

A 4 = 0.14 
A s = 0.12 


optimal control inputs are shown in 


Fig. 3.15 for four different values of k = Arj = fc 2 . The residual vibrations have 
not increased and the response is robust to the switching time. Now consider 
truncating the solut ion to one decimal place as: 

to = 0.0, Ao — 0.1 
tj = 6.6, Ai=0.1 

t 2 — 2.3, A 2 0.2 

t 3 = 9.0, A 3 = 0.2 (3.65) 

U = 4.7, A< = 0.1 
t$ = 11.4, As = 0.1 
t f = 11.5 

The time responses of x 3 to the robust time-optimal control inputs are shown in 
Fig. 3.16 for four different values of k = fcj = fc 2 . Again the residual vibrations 
have not increased but the overall rigid body motion has an offset. 


3.4 Conclusions 

Robust fuel- and time-optimal control problem of flexible spacecraft 
in the face of modeling uncertainty has been investigated. In particular, the 
rest-to-rest maneuvering problem of a flexible spacecraft equipped with on-off 
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reaction jets is considered. The rigid body solution is shown to have signifi- 
cant residual vibrations even for cases with prolonged maneuver time and it is 
necessary to include the flexible modes in the model. The results indicate that 
the proposed approach significantly reduces the residual structural vibrations 
caused by modeling uncertainty. 
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0£u*<l and -l£u*$0 

Figure 3.1: Case 3-1 



Figure 3.2: Case 3-2 
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Figure 3.8: Responses to fuel and time-optimal control input (Case 2-1). 
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Figure 3.14: Responses to fuel and time-optimal control input (Rigid Body 
Solution). 
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Figure 3.15: Responses to robust fuel and time-optimal control input (Case 







Chapter 4 


Sliding Mode Control Design for Uncertain Flexible 

Structures 


4.1 Introduction 

Open-loop controllers do not have feedback information of errors 
which arise due to uncertainties in the plant. As a result it would be bet- 
ter to use closed-loop feedback controller which should be robust to allow for 
plant uncertainties. Any feedback design which provides better robustness to 
parameter variations with no residual vibrations would be a viable alternative 
to robustified time optimal control logic. 

The ultimate goal of the future is to develop feedback control logic for 
achieving the robust time-optimal performance presented in the earlier chapter. 
Any scheme should not cause structural excitation and should be robust to 
parameter variations. The sliding mode control method is used in an attempt 
to design a controller with the desired properties. The idealized sliding mode 
control has a high frequency control activity and this is eliminated by the use of 
boundary layer. The sliding mode control allows for parameter variation of the 
nominal plant parameters. Since the sliding mode control method does not use 
bounded controls, a saturator is used to attain bounded controls. The boundary 
layer concept eliminates chattering by trading off with tracking performance. 
The good choice of the sliding function which is a weighted sum of the errors 
enables one to derive a control law. The sliding control method is applicable to 
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simple systems with decoupled equations of motion. The equations of motion of 
flexible spacecraft is modelled with the spring mass configuration. The modal 
equations are decoupled and this makes it possible to obtain the sliding mode 
controller. 

The remainder of this chapter is organized as follows. Section 3.2 
describes the variable structure system from which the sliding mode control is 
developed. Section 3.3 discusses the mathematical details of the sliding mode 
control theory. The boundary layer method and the robustness to parameter 
variations as applied to ideal sliding mode control is discussed. In Section 3.4 
the benchmark problem of an uncertain flexible structure is defined and the 
sliding mode control design is shown in Section 3.5. The simulations of the 
responses are shown for the nominal and perturbed parameters of the model. 


4.2 Variable Structure Control 


The variable structure control uses switching logic for control. This 
utilizes functions which would be used as switching functions for switching the 
control actions to on and off states. This switching function uses information 
of output feedback. 

Consider a system of the form : 


i = f{x,t) + B{x,t)u 


(4.1) 


where / and B are matrices and u is the control vector. The components of 
this vector u is determined by 


, / ufl *.(*) >0 (a o) 

“ i(l '‘ ) = U-i ».W<0 1 ] 

where s, ( i ) is the switching function and uf ,u t are the components of control 
u. The phase plane of the dynamics are utilized in designing the switching 
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functions to yield the control laws. The sliding mode control has been devel- 
oped from variable structure control. The variable structure control system 
has several substructures which may be unstable. Such cases are unacceptable. 
Sliding mode controller can be obtained as a special case of variable struc- 
ture control when the right hand side of the equations of motions describing 
the dynamics has discontinuities on certain surfaces in system phase spaces. 
The discontinuous control actions of the sliding mode can be viewed as an at- 
traction of phase trajectories onto the sliding surface in phase space. When 
the trajectory approaches the sliding surface, the perturbations about the sur- 
face becomes discontinuous causing high frequency chattering on/off of control 
states. This is highly undesirable. Sliding mode control provides robustness 
to parameter variations, robust performance to take care of plant uncertainties 
and external disturbances. 

4.3 Sliding Mode Control Theory 

The sliding mode control system is a high gain system and utilizes a 
high control authority for the tracking problem. The large bandwidth provides 
a fast settling time. Robust control using a feedback scheme is more effective in 
treating the plant parameter uncertainties as compared to open loop schemes. 
The sliding mode control is employed for designing the control law for uncertain 
flexible structures. The sliding mode control has been developed from the vari- 
able structure mode control [22,23,24,25]. The sliding mode control is robust 
to plant uncertainties and external disturbances. The trajectory of motion of 
system as illustrated in Fig. 3.1 is causing discontinuity in control action. This 
is the case when the trajectories are attracted from both regions Q + and Q~ . 
Once it reaches this discontinuous surface the motion occurs along the surface. 
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This motion is known as the ‘Sliding Mode’. This is the ideal case of sliding 
mode. The real sliding mode occurs with uncertain plants having parameter 
variations and unmodelled dynamics. The real sliding causes finite frequency 
control switches and the trajectory to jump across the discontinuous surface 
s,(x). Consider the system 

i = (4.3) 




( 


u, 

u” 


s,(x) > 0 
s,(x) < 0 


(4.4) 


where x = (ij, x n ] T are the coordinates of the system, / = [/i,-..,/n] T 

and u = [ui, . . . ,u m \ T . has discontinuities on the surface s,(x) = 0. 


A sliding inode exists for a region of non-zero measure on the surface 
s(x) = 0 if the projections of / + = /(x,t,u + ) and /“ = f(x,t,u~) are of 
opposite signs on the tangent to the surface and are directed towards it. Abo, 
we have 


lim i 

j—o- 

< 0 

(4.5) 

lim i 
»— 0+ 

> 0 

(4.6) 


The solutions of the system of equations with the right hand side discontinuous 
appears in the sliding mode. Filippov’s continuation method is used in this case. 
The system under consideration is linear with respect to u. 

Let G be m x n matrix whose rows are given by 


d . , d 



(*)j 


(4.7) 
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The equivalent coni rol u is obtained as 


Gf + GBu = 0 

u = -[GB}~'Gf 
x = f-B[GB]~ l Gf 


(4.8) 

(4.9) 
(4.10) 


provided [GB\~ l exists and the initial conditions s(x(0)) = 0. This is based on 
the following property 

i(x) = 0 (4.11) 


Now consider the rase of real sliding mode where the motion occurs near a 
finite 6 vicinity of the manifold s(x) = 0; 


II •{*) 


«<M -(*) 11= 


m 



(4.12) 


Let r(s,x) be the distance between any point in the neighbourhood and the 
manifold s = 0. It is given by 


r(s,x) < PS 

where P is a positive scalar. The real sliding case results in 


(4.13) 


x = f(x,t) + B(x,t)u(x,t) 


(4.14) 


where ti incorporates non- idealities and Bit satisfies the Lipschitz conditions 
(piecewise continuous and bounded). This is in agreement with most real 


systems. 


i = /°(*,<) 

f = m / + + (!-/•)/' 


(4.15) 

(4.16) 
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for 0 < /i < 1 and f° is a convex hull of f + and /“. Define f 0 , /+, /- as vectors 
joining origin to points /°,/ + ,/' respectively and gradient vector Vs. From 
Fig. 3.2 it can be s<«n that 

Vs-fo = 0 (4.17) 


Due to frequent switching let the state velocity be / + and f~ at time St i and 
St 2 respectively. Then substituting 


6t j 

^ St i •+• Sti 


(4.18) 


into Eqs. (4.16) we get 


f° = (SUf- +6t 2 f-)/(6t l +st 2 ) • (4.19) 


i.e., f° is the average velocity of motion. Now expanding Eqs. (4.17) we get 


[W + MI -fi)r}(d/dx)s = o 

11 (/- -r)(d/dx)s 

Substituting for /i in Eqs. (4.16) we get 


(4.20) 

(4.21) 


i = [(f-(d/dx)s)f + - (f + (d/dx)s)f-)/[(f- ~ f + )(d/dx)s) (4.22) 


This can be regarded as the sliding mode equations. 

Theorem 

If 


1. On the interval [0, T] any solution x(t) of the system described by 
Eq.(4.14) is such that the state trajectory is in 6 neighbourhood of the 
manifold s(x) = 0, or the inequality (4.12) is valid; 
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2. For the right hand side of Eq.(4.10) obtained using the equivalent control 
method 

X- = f{x\t) - *) ( 4 - 23 ) 

a Lipschitz constant / exist and x* is a solution of Eq.(4.14); 

3. Partial derivatives of the functions B(GB)~ l exist and are bounded in 
any bounded domain. 

4. For the right hand side of Eq.(4.14) there exist positive constants q and 
r such that 

II /( x, 0 + B(x , f)u(x, t) || < q + r || x || (4.24) 

then for any pair of solutions to Eq.(4.14) and Eq.(4.23) under the initial con- 
ditions 

|| x(0) — x*(0) ||< pS (4.25) 

there exists a positive quantity h such that 

|| *(0 -*“(<) ||< hS (4.26) 

for t € [0, T\. 

The proof the above theorem lies in showing that the ideal sliding 
and real sliding are close together provided the conditions stated hold. For a 
real sliding i ^ 0 as is the case in ideal sliding. The control u will be obtained 
as 

ti = -( GB)- l Gf + ( GB)~ l s (4.27) 

and substituting ti into Eq.(4.14) we get 

x — f(x,t) — B{x,t)[G(x)B(x,t)]~ l G(z)f(x,t) + B(Xyt)[G(x)B(x,t)] s 

(4.28) 
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The solution of Eq. ( 4.10) is given by 

*•(0 = x-(0) + jf {/(*•, r) - B( r)[C(*-)B(x-, r)]-‘G(x')/(x-, r)}<fr 

(4.29) 


The solution of Eq.(4.28) is given by 


x(t) = x(0)+ /‘{/(x,r)-B(i,r)[G(x)5(x,r)]- l G(i)/(x,r)}^4.30) 

+ / B(x, r)[G(x)B(x, r )] -1 idr (4.31) 

Jo 

Now subtracting Eq.(4.31) from Eq.(4.29) gives 

II x(0 -*•(<> II < pi+ /' Ml * - *■ II *■ + 

JO 

|B(x,T)[c(x)a(x,r)|-'* in; 

+ 1' || ±B{z,t)IG(z)B{x,t))-' nil 4 II dr (4.32) 

The condition given by Eq.(4.24) and the bounded solution of Eq.(4.31) can 
be combined as 

II x«) ||<|| x(0) II +mr+ f T n II x || dt (4.33) 

JO 

using the Bellman-Gronwall lemma, the above inequality is 


II *(<) ll< (II 1(0) || +m7>" r (4.34) 


By using the conditions 1-3 of the theorem and inequality Eq.(4.34), we get 

|| x(t) - x*(t) ||< p£ + l f || x - x' || dr (4.35) 

Jo 

where p is a positive scalar. Bellman - Gronwall lemma applied on the above 
inequality gives 


l|x(t)-x-(l)||<M 


(4.36) 
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where ( t 6 [0, T]) and h = pe IT . This shows that if initially the actual and 
ideal states are sufficiently close together, then they remain close over a finite 
time. This result could also be extended to infinite time interval. 

Now consider the exponential stability of linear time invariant systems 
described by 

i = Ax + Bu + g (4.37) 

where A and B are constant matrices and disturbance is represented by g. The 
sliding mode occurs along the intersections of s, = 0 and satisfies 

^ = CAx + CBu + Cg = 0 (4.38) 

This results in equivalent control as 

* = -{CB)~ x {CAx + Cg) (4.39) 

which becomes 

X* r: [/ - B{CB)~ l C)Ax m + [/ - B{CB)-'C\g (4.40) 

where x* is the variable of ideal sliding case. In the case of real sliding case the 
above relations are 

ti = -(C'SrHAx + ^ + tCB)'^ (4.41) 

x = (/- B(CB)-'C]Ax + [I- B{CB)-'C\g +B{CB)-'s (4.42) 

The solution of the above two cases can be obtained as : 

x*(<) = *(t)x'{0)- J*9(t-r)[I -B{CB)- l C)g(r)dT (4.43) 

x(t) = $(t)x(0) — /* $(t — t)[I — B(CB)~ l C]g(r)dT 

Jo 

+ j\(t-T)B{CB)- x j^dT 


(4.44) 
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where 4> is the state transition matrix. Integrating the above and subtracting 
them and taking the norms of both sides gives 

II x(0 - *•(£) II < II m || || x(0) - x*(0) II 4- II B || || (C5)- 1 II || s(t) || 

+ II B || || (CB)~ l || || *(<) || || S (0) || 

+ £ II jr*(( - r)B(CB)-' || || * || dr (4.45) 

For real sliding the inequality || s ||< 8 holds and the functions 
|| $(<) || and /o || j;$(t — t) || dr are bounded in the infinite time inter- 

val for asymptotically stable systems. Then the following inequality holds for 
positive n 

||x(*)-x*(*)||<n* (4.46) 

and 

limx(t) = x*(f) (4.47) 

for 0 < t < oo. 

Now the F ilippov’s method is applied to the above problem of linear 
time invariant systems. The convex hull in this case is given by 

+ <)£>,“, (4-48) 

1=1 

where /i, > 0 and Pi = 1* The velocity vector lies on the sliding surface 
hull and hence Gf° = 0. Let z — £*7, P> u i and substituting it gives 

G{f + Bz) = 0 (4.49) 


or 


2 = -(GB)~ l Gf 


(4.50) 
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which is the identical to Eq.(4.10). Thus the Filipov’s method is in agreement 
to the equivalent control method for the linear time invariant systems. The 
following conditions hold for attraction to the sliding surface 


f + (d/dx)s < 0 ,x€S + 

(4.51) 

f~(d/dx)s > 0 ,x€<7” 

(4.52) 

which can be combined as 


j t $ 2 {x\t) < 0 

. (4.53) 

which is a local sliding condition. For a proper function s(x,f) and a positive 

scalar rj we have 


^s 2 (x;t) < -r}\s\ 

(4.54) 

as the global sliding condition for a time varying case. Consider the 
single-output linear time-varying system 

single- input 

x[ n) + a„_i(t)x ( 1 n_1) + . . . + a 0 (t)*i = « 

(4.55) 

Now define the states as x(t) and the desired states xjj, ) as 


x(t) = [x,(f),x 1 ,...,x ( i n-,) (t)] 

(4.56) 

x d {t) = [xii(t),x rfl ,...,xjj _1) (t)] 

(4.57) 

X = x(t) — Xd(t) 

(4.58) 


Assume 

1. x(0) = 0 

2. |xi n+,, (f)| < r for all time and constant v 
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Define the time varying S{x;t), the sliding surface as 

s(x;t) = Cx(t) = 0 (4.59) 

where C = [c,, (^. 1 , 1 ]. Now the objective is to obtain control u(t), to 

make s{x\t) = 0. Expanding the above equation we get 

+ "f = xir" + £<*,*1? (4.60) 

1=0 1=0 

Assumption 1 and t he uniqueness of solution of ordinary differential equations 
would result in 

x(t) = x d (t)Vt (4.61) 

Set u as 

a = 0 T (x)x + ki(x\ t)ii + 1 - k n sgn{s) (4.62) 

1=1 

Substituting for control in Eq.(4.55) to get 

“.s 2 (x; t) ~ £(£,(x) - a,-i(f))x,s + ^(c, + *r,(x; t))x, +1 s 
» = 1 

-sx d<n+1 - fc„|s| (4.63) 


The sliding condition results in the constraints 

0i{x) = < ai.i(t),XiS > 0 V*,i = (4.64) 

0i{x) = 0~ > ai.i(t) y XiS < 0 Vi,* = l,...,n (4.65) 

ki(x]t) = k+ <— Ci,x, + i3>0 Vt,*‘ = 1, . . . ,n — 1 (4.66) 

ki{x\t) = k~ > c,,x <+ |5 < 0 Vt,i = l,...,n - 1 (4.67) 

K > v (4-68) 

When c, = 0 for some i, corresponding fcjXj+i is set to zero. The discontinuities 
are at 

x, = 0; ij = 0 ,s = Cx = 0 (4.69) 
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but control is discontinuous only at s(t) = 0. Combining Eq(4.54) and Eq(4.63) 

—»*(*! <)<-(*.- »)l»l (4-70) 

for all t. 

Parameter variations in sliding mode control is incorporated in the 
constraints. Consider the finite bounds of parameter variations of a, as 


o, < <*,(0 < 7ii t = 0, . . . , n (4.71) 

Then from Eq(4.64)-Eq(4.67) we get 

ft < *.-» (4-72) 

ft > 7«-i (4.73) 

for i = l,...,n, control yields x(t) = xj(f). To incorporate robustness to 
external disturbance given by 

d(x\t) = [0, . . . , 0, di[x,t)] T (4.74) 

\di(x\ <)| < £«.W + <o (4.75) 

1=1 

The control law given by Eq(4.63) will yield x(t) = x,i(t) as long as 

ft < or,-i - Si (4.76) 

ft > 7<_i -6„t = (4.77) 

k n > v + 6 0 (4.78) 


are satisfied. The initial condition assumption, xj(0) ^ x(0), may be true and 
this implies that the points in hyperspace may not be on the sliding surface. 
But using the control law, the offset will asymptotically approach zero provided 
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Eq(4.60) is stable. Let the class of nonlinear time-varying systems described 

by 

^ T ‘ J ' = / ; (0i,0 2 ,...,0 p ;O + j = l,...,p (4.79) 

where, for : = 1,. . . ,p and 0, = [0,, 0,, . . . The tracking problem is 

to make 6j approach 6^. Define the tracking error as 0,(f) = 0,(t) — 0<<,(O- 
Now replace the discontinuous control by the continuous control law 

5,(0,; t)mCA(t) (4.80) 

Define a boundary layer about the sliding surface as 



(4.81) 


(4.82) 

where C Jt > 0 and e, is the hyper-radius of boundary layer, 
layer Bj(t) is defined by 

The boundary 

Bj(t ) = {© : sJ(Q i; t) > 0| jt(0 i; O < 0} 

(4.83) 

d _ d d . 

dt 3j ~ dt S} ~ dt Sj 

Choose Uj(Q,t) so that outside B } we have 

(4.84) 

j t >; > ove € s;(t) 

(4.85) 

j t >* < ove €«+(<) 

(4.86) 


The boundary layer is the attraction for points in hyperspace which is assured 
by the above conditions. Any continuous interpolations between u, for s~ and 
Uj for s + will suffice. If 


(4.87) 
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then with 0^(0) = 0<*,(O) 

|©,- -0*| < tjVt >0 (4.88) 

If ©j(0) ^ 04 ,( 0 ), we have 

|Bj - ©*| < + P(t) || 0,(0) || exp(-At) (4.89) 

where P(t ) is a polynomial in t. 

Consider t he single-input single-output system of the form 


i, = x 2 

it = *3 

in = /(*) + $(x)u + d (4.90) 

where the states are i = [ij z 2 ... x n ] r , u is the control, y is the output and d 
the disturbance. The tracking problem is to make output approach the desired 
output while allowing for uncertainties in / and g. The time varying sliding 
surface s is defined as 

*(*;«) = (^ + A ) B-,i > A >° ( 4 - 91 ) 

where e = — y and A is constant scalar. The trajectories are attracted to 

the sliding surface and this results in nullifying the error. The phase plane 
trajectory should converge to the sliding surface in order for the control law to 
be effective in fulfilling the desired tracking objective. Chattering is removed 
by smoothing out ihe control surface discontinuity in a thin boundary layer 
neighbouring the switching surface which is given by 

B(t) = {x : |«(x;l)| < A""'e} 


(4.92) 
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where the thickness of the boundary layer c > 0. Consider the modelling error 


/ = / + */ 


(4.93) 


Assume that 


\6f\ < F (4.94) 

|d| < D (4.95) 

where /, F and D are known functions of time and states. Let 

k(x; 1 ) = h(x,t)+d(x,t) + v(t) + rj >r) >0 (4.96) 

a = -f -Y^Cp~ l X p x n ~ p - k{x-,t)sgn(x) (4.97) 

p=i 

Introduce a boundary layer of thickness e, and thereby modifying the control 
to 

u = -f - £ C;- 1 \ p x n ~ p - k sat(j^ -) (4.98) 

p=i 

where sat is the saturator and this control does satisfy the sliding condition 
Eqs.( 4.54). For points satisfying x(0) = 0, 

i = M*i Ojs; + (</(*■ O + ^O-*?’) ( 4 -"> 

= + + + (4.100) 


for 6 f and k are continuous functions. 

The sliding condition that needs to be satisfied for convergence is 

< -»MI (4101) 


for a non-negative constant r j. This condition is like a Lyapunov function. It 
shows that the distance squared from the sliding surface ( s 2 term), gets smaller 
as the trajectories are traversed toward the sliding surface. The sliding surface 
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also allows for disturbances and uncertainties like parameter variations and 
the sliding condition takes care of the perturbation caused. Control laws which 
satisfy the sliding condition are discontinuous across the sliding surface and 
lead to control chattering. This is due to the absence of the exact knowledge of 
the value of the sliding variable. This leads to jumps across the sliding surface 
at high frequency. This is the phenomenon behind high frequency chattering. 
Chattering can excite higher order modes. Thus it lowers the robustness to 
unmodelled dynamics which are the higher modes. Hence it would be desirable 
to smoothen the control and eliminate chattering. 

In case of g{x) = b (scalar constant), let b vary in the interval given 
by 

0 < b mtn < b < b max (4.102) 

This bounded interval margin is assumed to contain the parameter 

variations of parameter b. This should take care of robustness to parameter 
variations in 6, the control gain constant that appears in the equations of mo- 
tion of the system. This interval margin is bounded and represents a function 
of masses and spring stiffness in the case of flexible structures. Hence the vari- 
ations in b are caused by variations of real parameter variations like masses 
and spring stiffness which are uncertain parameters. The control gain b ap- 
pears as a product in the equations of motion and this results in multiplicative 
uncertainty. Define 

6 = {bminbmax)* (4-103) 

0 = (4-104) 

is the analogue of gain margin in linear control and 6 is the geometric mean 
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of extreme bounds of b. The perfection in tracking is achieved by 

A n e = 0 k ma£ (4.105) 

In the case of time varying boundary layer the control is chosen to satisfy 

+ A-*|m|W - 1) (4.106) 

2 (it 

The corresponding dynamics of boundary layer thickness are 

k(x d -,1)>\ n t/0 => e + Xt = jh k (x*\t) (4-107) 

k{x d ;l) < \ n t/0 => c + — = 0 (4.108) 

e (0) = ^k(x d ( 0),0) (4.109) 

k(x;t) = [k(x-,t)-k(x d ;t)] + \ n e(t)/l3 (4.110) 


The control law is given by 

u = h~ l [u — ksgn(s)] (4.111) 

The sliding condition is satisfied and this represents the analogue of the switch- 
ing function of variable structure control. The control law takes into account 
parameter variations in b. Chattering is greater and it would be desirable to 
smoothen the control action. This smoothening should keep useful properties 
of robustness to parameter variations, robustness performance and good track- 
ing properties of sliding mode control. Slotine and Sastry [22] have developed 
a boundary layer to eliminate chattering. Redefining the boundary layer <j> as 

<p(t) = A n_1 (t)c(t) (4.112) 


modifies the control to the following form 


u 


-f - £ C; -1 A ,, - 1 (Ax (B - ,,) + pAx n-,>-1 ) 

pml 

-T(x;t) sot(^) 

C-b 


(4.113) 
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The boundary layer is space surrounding the sliding surface in am envelope 
in the hyperspace. This contains the discontinuity jumps which take place 
and could be thought of as an analogue of a dead zone. The undesirable 
chattering is eliminated. The action of boundary layer is similar to that of 
a low pass filter near the sliding surface to variable s. A trade off between 
tracking and robustness to parameter variations results. The tuning parameter 
is the thickness of the boundary layer. A time- varying boundary layer is utilized 
to allow for variations of boundary layer with time. In the case of time varying 
boundary layer, the dynamics are given by 

k(x d ) > \<j>/0 d =» 4> + X <f> = 0 d k(x d ) (4.114) 

Lix d ) < * j> + \<f>/fi = k(x d )/0 d (4.115) 

where 

T as k(x) - k(x d ) + A <t>ff3 d (4.116) 

and x d is the desired target state x. The initial condition is 

*(0) = 0 d k(x d (O))/\ (4.117) 


4.4 Model of an Uncertain Flexible Structure 


Consider ihe two mass-spring system shown in Fig. 3.3 which is a 
generic model of an uncertain dynamical system with one rigid body mode and 
one vibration mode [26] . The equation of motions can be written as : 


■ Xi ' 


■ 0 0 10' 


' x x ' 

±2 


0 0 0 1 


Xj 

*3 


—k/mi Jfc/mj 0 0 


*3 

. *4 . 


fc/m 2 —kjmi 0 0 . 


. x* . 


0 

0 

l/ m ! 

0 


u (4.118) 


y = *2 


(4.119) 
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where x x and x 2 are the positions of body 1 and body 2, respectively; x 3 and 
x 4 are the velocity of body 1 and body 2 respectively. The aim is to design a 
feedback controller for a unit step output command tracking problem for the 
controlled output x 2 with the following specifications : 

1. Control input is limited as |u| < 1. 

2. Performance requirements : settling time and overshoot are both to be min- 
imized. 

3. Performance robustness and stability robustness with respect to the uncer- 
tain parameters m, . m 2 and k ( with nominal values of mj = m 2 = k = 1) are 
to be maximized. 


The equations of motion can be written in modal coordinates as : 


y, = <f> x u 
h + w 3 y 2 = 4>2 u 


(4.120) 

(4.121) 


where u> is the natural frequency of the vibration mode and <f > i , are the gain 
constants. 


4.5 Sliding Mode Controller Design 

The initial conditions for the rest-to-rest maneuver is 

yi(0) = y 2 (0) = 0 
y'i(0) = y'j(0) = 0 

and the final conditions 

yid = 1) y 2 <i — 0 

ilid = o, y u = 0 


(4.122) 

(4.123) 


(4.124) 

(4.125) 



90 


The sliding surface is to be defined in terms of sliding variable s. Let 

s = y x + A y 2 (4.126) 

where y = y — yd- The parameter A is the weight used in the error function 
and it also acts as the poles of the filter in the dynamics of the sliding variable. 
Thus the problem is reduced to a tracking problem to satisfy the rest-to-rest 
maneuver. By approaching the sliding surface in the phase space would be the 
strategy of the coni roller for achieving perfect tracking. Differentiating s, we 
get 

,s = — Aw 2 y 2 — jjij — Ay 2 j (4.127) 

where yd is the desired target acceleration. 

The control logic u is of the form 

u = a{u — fcs<7n(s)} (4.128) 

where 

u = a[Au> 2 y 2 + Ay 2 << + yu] (4.129) 

a = (^i + A^) 1 (4.130) 

and the Cj,4>\,4>i are the worst case parameter variations. This is introduced 
to make the control robust to parameter variations in masses ml and m2 and 
spring stiffness k. 

To satisfy the sliding mode condition 

si < — »j|s| (4.131) 

the following condit ion need to hold 

k > j— [— // + |Au/ 2 y 2 + yu + Ay 2l *|] + |aAu> 2 y 2 | + |A(yi,* Ay 2 j)| (4.132) 

<f>i + A fa 



n 
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This sets the lower bound on the value of K to be chosen. Now the boundary 
layer 4> is included io modify the control law to 

u = {u ~^sat(s/4>)} (4.133) 

To satisfy the bounded control specification we saturate it as 

u *- sat(u) (4.134) 

Carrying out the simulations the results for the above example yields a sluggish 
response. The response has considerable fluctuations and would result in large 
settling time larger than 10 seconds. It shows some chattering activity initially 
which is undesirable. By using 

5 = yi + Ay 2 (4.135) 

and using the same control law obtained before we obtain the better response 
shown in Fig. 3.4. This design has settling time of about 8 seconds and is 
simulated for the exact nominal case of parameters. It has a small overshoot 
and it has no chattering or any residual vibrations. The control time history 
also shows the switching of control to both extremes and also the linear varia- 
tions combinations. The rise time is fast and these are the desirable properties 
required for a control design. 

This control design should be robust to parameter variations. Now 
the parameter variations are introduced and the simulation results are shown 
for various values of parameter variations in the Fig. 3.5-3.8. This illustrates 
the robustness to parameter variation and stability robustness. The ability of 
the controller to switch and provide the necessary control action is seen clearly 
in these simulations. Further the variations in the control actions to different 
cases of parameter variations are seen. This is in contrast to the robustified 
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time optimal solutions where one control time history had to be used for any 
case of parameter variations. 

4.6 Conclusion and Issues 

An observor can be designed for estimating one of the states the 
position of the first mass. However the methodology used in the design of the 
sliding mode controller is to convert to modal coordinates. The transformation 
from read coordinates to modal coordinates needs the knowledge of both the 
states. The sliding mode controller is a nonlinear scheme which also provides 
robustness to parameter variations and robust performance. Its introduction 
to a generic case of uncertain flexible structures has been shown to provide a 
good performance. 
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Figure 4.4: Response and control of sliding mode controller for the case 
with parameters k =1 , ml = 1 , m2 = 1 (nominal case). 
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Chapter 5 


Conclusion 


The robust time-optimal open-loop control problem of flexible space- 
craft in the face of modeling uncertainty has been investigated. The primary 
study objective was to explore the feasibility of computing open-loop, on-off 
pulse control logic for uncertain flexible spacecraft. The results indicate that 
the proposed approach significantly reduces the residual structural vibrations 
caused by modeling uncertainty. The results also indicate the importance of a 
proper jet placement for practical trade-offs among the maneuvering time, fuel 
consumption, and performance robustness. It is shown that most problems re- 
quire an accurate mathematical model and thus are sensitive model parameter 
variations. A nonlinear constrained parameter optimization method is used for 
solving this problem. The structural flexibility and mass distribution of the 
vehicle is uncertain while the total mass is well known. So the uncertainty 
occurs in modal frequency and mode shapes. The case of one sided control 
inputs only satisfies the necessary conditions and theorems on uniqueness and 
structure of controls do not exist at this time. The time-optimal responses are 
shown to be quite sensitive to variations in model parameters. The response of 
displacement of masses due to variation of spring stiffness is illustrated. Sim- 
ilar responses can also be observed for arbitrarily combined variations of k, 
and m, but keeping total mass constant. The two one-sided control inputs 
solutions are bang-off- bang nature and in this case the control on-time is not 
equal to the maneuver time as was the case in bang-bang control. This shows 
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saving in fuel consumption even though the objective function does not have 
any weighted fuel terms. The robustness constraints are incorporated and the 
parameter optimization problem is solved with these constraints. The number 
of pulses are increased to match the additional number of constraints. The 
results indicate that the robust time-optimal solution are less sensitive to pa- 
rameter variations. The performance robustness increases at the expense of 
increased maneuver time. By prolonging maneuver time does not reduce resid- 
ual vibrations caus^rd by model uncertainty but a properly coordinated pulse 
sequences, the flexible modes are not significantly excited during maneuver and 
are well desensitized. The summary in chapter 2 shows that case 2-3 was the 
best actuator location in terms of overall performance. The accuracy of the 
switching time is also presented and it is shown that the robust time-optimal 
solution is insensitive to parameter variations for switching time accuracy upto 
two places of decimal. The problem of solving the optimal actuator location 
using m thrusters in an n mass spring configuration is still an open research 
problem. 

The robust fuel- and time-optimal control of uncertain flexible space- 
craft have also been investigated. The consumed fuel was assumed to be small 
compared to the mass of the spacecraft and it was assumed to be proportional 
to the jet on-time. The performance index is a weighted function of the ma- 
neuver time and jet on-time. The bang-off-bang nature of the two one sided 
control provides a considerable reduction of jet on-time. It is again noticed 
in the fuel- and time-optimal control solution that the total maneuver time in 
not the same as the jet on-time. The rigid body solution is shown to have sig- 
nificant residual vibration and hence it is necessary to include flexible modes. 
The robust fuel- and time-optimal control again shows increased maneuvering 
time at the expense of robustness. It was shown that by obtaining a rigid body 
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solution with comparable final time with the robust fuel- and time-optimal so- 
lution of flexible structure, the residual vibrations are significant larger in the 
rigid body solutions. 

The sliding mode control design using bounded control inputs was 
obtained for the benchmark problem. The results indicate robust performance 
to parameter variations. The feedback design uses a saturator to bound the 
control inputs. The stability of such controller was not studied and it remains 
an open issue. This scheme is not minimizing the maneuvering time or the fuel 
consumption. The fast settling time and low overshoot are the characterstics 
shown by this design. 

There are some limitations of the present work. The analysis is re- 
stricted to linear, elastic and undamped systems. The future research directions 
should be towards i he ultimate goal to develop nonlinear feedback control for 
achieving the robust fuel- and time-optimal performance. The experimental 
verification of the applicability of the robust fuel- and time-optimal control 
should provide more insight in its practical implementation. The mathemati- 
cal rigorous conditions are also needed for necessary and sufficient conditions 
of time-optimal control in the case of two one-sided controls. 
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Appendix A 


Numerical Solution of Nonlinear Constrained 
Optimization Problem 


The constrained parameter optimization problem can be solved using 
numerical schemes [1 8]- [21}- The problem can be stated as 


mtn /(x), x € R n 


9j( x ) = 0, ; 
9j( x ) > 0, j 


1 , . . . , m e 
m e + 1, . . . , m 


X; < x < x„ 


(A.l) 


where x/ and x tt are the lower and upper bounds of x respectively. It is also 
assumed that the functions are all continuously differentiable. 

Let the current iterate be x* and v k be the corresponding approxima- 
tion of optimal Lagrangian multiplier. The positive definite approximation of 
the Hessian, B k , of the Lagrangian function 

L(x, u) = /(x) - £ u i9j( x ) ( A-2) 

j=i 

where u = (ui, € R m> ■ Let the value m' = m + 2n. Let the bounds 

of x be defined as 


gj{x) = z b" m ) _ j = m + 1, . • ■ , m + n 

g } (x) ---- x[r m - n) - j = m + n + 1 m' (A.3) 


Linearizing the constraints of the optimization and minimizing quadratic ap- 
proximation of the Lagrangian function as 


min ^ d T B k d + Vf(x k ) T d (A. 4) 

Vtf j(x k ) T d + g } (x k ) = Q,j = 

Vfjj(x k ) T d + gj{x k ) > 0,j = m e + 
x/ — x k < d < x n — Xfc 

Let d k be the solution and u* be the corresponding vector of the Lagrangian 
multiplier. 

x k +i = x k + a k d k (A. 5) 

where a* is the line* search parameter such that it decreases the ‘merit function’. 
However Eq. (A.4) may lead into a non-feasible problem but it may still be 
feasible due to linearization of the constraints. Also it needs the gradients 
of the constraints at each iteration. To avoid these problems Schittkowski 
proposes an augmented Lagrangian method using the an additional variable 6, 

mm ^ d T B k d + f(x k ) T d + ^pkS 2 (A.6) 

Vgj(z k ) T d + (1 - %j(z*) > OJ € J k . 

?9j(*kti)) T d + gj(x k ) > 0 ,j € K k 
X l < X < x u 

0 < S < 1 , 


where 


J k = { 1 , . . . , m e } U 

{> :m t < j < m,gj(x k ) < t or > 0}, 


(A.7) 


Ill 


where e is the accuracy required and pk is the additonal parameter used to 
reduce the influence of 6 on the solution. The matrix B k update can be used 
from the BFGS update formula given by 


Bk+x 

Pk 

<Jk 


R > qkq k _ BkPkPkBk 
k q£qk Pk B kPk 
x k + 1 %k 


v/(n +1 ) r - V/( z k ) T 


(A.9) 


The constraints are non-linear and the cost function is linear. The 
standard IMSL subroutine NCONF uses successive quadratic programming 
method to solve a general nonlinear programming problem. It requires the 
value of the constraints and their gradients at each iteration. This scheme 
uses an initial guess value and then iterates to converge to the solution. There 
are many local optimum solutions which are obtained by using different initial 
guesses. However the optimum solution is one which has optimal cost function 
amongst the local optimum solutions. 




